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ABSTRACT 

The equivalent width dis tribution f unctio n (EWDF) of H I absorbers specific to the void environment 
has been recently derived ((Manning! . l2002ft . The findings revealed void line densities dM jdz ~ 500 at 
equivalent widths W > 15.8 mA (Nhi ;> 2.6 x 10 12 cm -2 ). I show that the void absorbers cannot be 
diffuse (or so-called filamentary) clouds, expanding with the Hubble flow, as suggested by N-body/hydro 
simulations. Absorbers are here modeled as the baryonic remnants of sub-galactic perturbations that 
have expanded away from their dark halos in response to reionization at z « 6.5. A 1-D Lagrangian 
hydro/gravity code is used to follow the dynamic evolution and ionization structure of the baryonic clouds 
for a range of halo circular velocities. The simulation products at z — can be combined according to 
various models of the halo vel ocity distributio n function to form a column density spectrum that can be 
compared with the observed ijMannind . l2002ft . To explain the observations with these models requires 
a search of parameter space somewha t beyond the envelope of conven t ion. F or a given circular velocity, 
a halo model more massive than the lNavarro. Frenk. fc White! 1199a, 119971 (NFW) halo is required to 
reproduce the observed line density of absorbers. A more massive, non-singular isothermal halo is used 
with a more favorable outcome. I find that such clouds may explain the observed EWDF if the halo 
velocity distribution function is as steep as that advanced bv iKlvpin et al . (1999). Observations are 
best explained when individual halos have sub-halos that occupy the flanks of the "parent" . A picture 
emerges in which growth by accretion of sub-halos is possible. Further analysis suggests that the mass 
distribution about a cloud may extend significantly farther than the virial radius. Accounting for the 
total void mass density remains an outstanding problem. 

Subject headings: intergalactic medium — quasars:absorption lines - dark matter - galaxies: halos 



1. INTRODUCTION 

Among the illuminati of cosmology, a paper that pro- 
poses to study Lya clouds in voids, and further seeks to 
model them as remnants of sub-galactic perturbations, is 
likely to be viewed as an anomaly; N-body/hydro simu- 
lations have had such succes s in reproducing the detailed 
nature of the Lya forest fe.p..lMir alda-Escude et ailll996t 
[ Hernauist et aUll996tlRle diger et al.lll99 8HCen fc Ostrikerl 
Il999t iDave et all Il999t iDave fc Trippl 12001ft . that seri- 
ously considering a different picture is nearly unthinkable. 

In these simulations, Lya absorbers are not physi- 
cally associated with discrete small halos; rather, the-v 
arise in large d iffuse sheets or filaments (e.g.. iBil 1199 3 
Weinberg etaLlll997tlMiralda-Escude et allll996HDave et all 
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McLin et all 12002ft. suggestir 



This view is sup ported by the conclusions of 




Dinshaw et al.l (|1997l ll998T) regarding the hit/miss statis- 



tics of moderate redshift Lya absorbers in double or group 
quasars. The possibility that these hit/miss statistics 
could instead be the signature of the clustering of small 
clouds w as considered, and reje cted without a thorough 
analysis ({Weinberg et all 11997ft 1 . Nor are Lya clouds 
expected to be found in voids at z w in any signifi- 
cant numbers, for in voids, the diffuse filaments expand 
with the Hubble flow and disperse to be low the current 
limits of detectability at redshifts z <; 1 l)Riediger et all 
1998J IDave et all 11999ft . Nevertheless, there exists now 
some considerable history of detections of low-redshift Lya 

1 However, in unpublished work, Manning (2000) produced a good 
fit to the hit/miss statistics over a range of 10 kpc to 1 Mpc using 
discrete Lya clouds clustered around small galaxy groups of physical 
extent ~ 750 kpc (e.g., "Local Groups"). 



are not as rare as the simulations would imply. Likewise, 
the discovery of compact h igh velocity clouds (CH VCs) in 
the galactic neighborhood <|Braun fc Burtorlll999ft . which 
are successfully m o deled as partially held by dark matter 
(|Blitz et all Il999t ISternberg et all 12002ft . suggests that 
sub-galactic halos with associated H I gas may not be as 
unlikely as previously thought. 

At issue is not just the specific disposition of dark mat- 
ter and baryons in vast, dark, under-dense regions. The 
observations of Lya clouds at low-redshift, where complete 
galaxy catalogs can be used to calculate their relative isola- 
tion, may provide a crucial test of the N-body simulations 
that play such a central role in informing our picture of 
process on the cosmic scale. 

Simulations have effectively been calibrated by the use 
of the high-redshift Lya forest. Under the assumption 
that dark matter (DM) can be char acterized by a 3- 
dimensional power spectrum P(k) (fsH Il993t ICroft et all 
11998ft . and that gas traces DM except at the smallest 
scal es, wh ere pressure forces are thought to be important 
|b1 Il993|) . P(k) is derived from high-redshift Lya forest 
spectra, and used to model the initial conditions of sim- 
ulations. The Lya forest thus became a testing ground 
for simulations, enabling researchers to gradually refine 
their code to a point that there is now excellent agreement 
with observations of a wid e range of Lya forest proper 
ties dMiralda-Escude et all 119961: iHernauist et all Tl 996: 



iRiediger et all 119981: ICen fc Ostrikeii Il999t IDave et al 



2 



Modeling the Void H I Column Density Spectrum 



fe Tripp| . l2001j) . These simulations show struc- 
ture formation and movement of matter from low, to high 
density region s in response to peculia r gravitational fields. 
For instance, ICen fc Ostrikerl l)1999f) use an open CDM 
model to follow the transfcrral of matter from underdense, 
to overdense regions, showing that over 90% of the mass 
in overdense (which we will later refer to as "shocked") 
regions at z = arrived after z — 3. 

The challenging task of following the movement of gas in 
a simulated cosmological context is handled by concentrat- 
ing baryonic and dark mass into baryon and dark matter 
"particles" , typically of mass ~ 10 8 , and ~ 10 9 Mg u re- 
spect ively ijKatz et allll996HHernauist et al.Lfl996tlDave et 
1999). Due to this coarseness, small galactic and sub- 
galactic halos cannot be represented. Thus, the baryons 
that could have been associated with smaller halos are 
distributed in a smooth manner through the large dif- 
fuse structures. However, the response of the particles 
to gravity produces convergent motions which are inter- 
preted as absorption lines - an effect known as t he "fluc- 
tuating Gunn-Peterson effect" l|Croft et allll99d) . Artifi- 
cial spectra are typically extracted by smoothing groups 
of particles with a spline k ernel with a radial scale that is 
dependent on the density ijKatz et all 11996(1 so that gas 
resolution varies from ~ 5 to J> 200 kpc for high, and 
low density regions, respectively l)Hernauist et all fl*996). 
Thus, even if sub-galactic halos could be resolved, the dy- 
namics of the behavior of baryons in the halo could not 
be realistically represented since the distribution of gas is 
determined completely by the spline function, which has 
no physical foundation. 

In attempting to characterize the nature of the evo- 
lution of structure seen in simulat ions, two broad types 
of environments are detected. In iRiediger et alJ Jl998) 
these environments ar e characterized as "sho cked" and 
"unshocked", while in ICen fc Ostrikerl l)1999(K they are 
"warm/hot", and "warm". In either case, the filling fac- 
tors, and the physics contributing to the physical state 
of the gas are approximately the same. The clouds in 
warm/hot, or shocked regions are collapsing structures, 
and are shock heated, while the warm, unshocked clouds 
characterize underdense regions in which gas is expanding 
with the Hubble flow, and is hence adiabatically cooled. In 
ICen fc Ostrikerl l)1999(l . the collapsing regions are thought 
to comprise about 10% of the volume of the universe at 
the present time. In the simulations of iDave et alJ |l999) 
and lDave fc Trippl l|2001j) . clouds at high redshift populate 
both overdense and underdense regions, but for z < 1, ab- 
sorbers are thought to arise in regions that are overdense 
by factors of 12 to 100; stronger lines originating close to 
galaxies, and weaker lines distributed in the diffuse fila- 
mentary structure. 

At z ~ 3, the situation is different; the apparent lack 
of a modulation in the incidence of Lya absorbers on ve- 
locity scales appropriate to galax y superclusters o r voids 
in high-redshift spectra prompted ICarswell fc Reesl l|1987|) 
to conclude that voids in the Lya cloud distribution can 
fill no more than 5% of the universe at (z) ~ 3.2, a result 
corroborated by oth e rs. H owever, the area is murky (see 
the review bv lRaucbl 1^998) for a discussion) since how one 
defines a Lya cloud void depends on an arbitrary detection 
threshold. 



However, the movement of this proto-void material to- 
ward mass concentratio ns will tend to e vacua te them. 
Simulations of voids by iMathis fc White! ((2002) suggest 
that no galaxy of any size is expected to currently popu- 
late voids of size ~ 10 h _1 Mpc. However, studies show 
that there are much larger voids, with an average diam- 
eter greater than twice as large (e.g., iHovle fc Vogelevl 
200 21 lEl-Ad fc Piranl l2000t). But void simulations of 
lArbabi-Bidgoli fc Mullen ((2002) appear to imply that larger 
voids have generally lower densities. However, the rela- 
tively large galax y density within the huge Bootes void 
(|Dev et all Il995|) , and the subsequent detection of other 
aQ, galaxies in 21 cm emission l|Szomoru sug- 
gests that larger voids may retain halos when smaller voids 
may not. 

The assumption that baryons trace a diffuse DM dis- 
tribution, and not smaller, sub-galactic potentials, can, 
in principle be tested by the search for Lya absorbers in 
local voids. For, if absorbers are diffuse, they will have 
expanded with the Hubble flow and now be so diffuse that 
they are unobservable (e.g. ?), while if they are held by 
sub-galactic halos, the gravitational potentials may be suf- 
ficient to restrain the evaporation of baryons and produce 
absorbers. 

Observations of low-redshift Lya clouds, aided by galaxy 
catalogs, can be used to determine the relative isolation of 
clouds, and may thus help to resolve these issues. For 
instance, Hubble Space Telescope (HST) spectra have de- 
tected numerous H I clouds for whic h no nearby galaxy 
was found llMorris et all Il991i Il993t IStocke et all Il995t 

81 iMcLin et all liF" 



iShull et all 119961: iTripp et all 119981 iMcLin et all I2002D . 

However, the methodology used by these investigators, 
though useful for gaining insight into the low redshift uni- 
verse, does not facilitate the study of void clouds in any 
systematic way. For, "nearest neighbor" measures of a 
cloud's isolation weights galaxies equally, whether dwarf 
or giant, and so lacks the sensitivity to mass that is intu- 
itively attached to the concept of isolation. 

A new method for resolving the e nvironments of Lya 
clouds is presented in lMannind l)2002|) (hereafter Paper 1). 
It uses ambient tidal fields, calculated from galaxy cata- 
logs, as a relative measure of a cloud's isolation, and thus 
weights the effects of galaxies by their mass, and by the 
inverse cube of the distance. Clouds can then be sorted 
by their degrees of isolation. A catalog of clouds can, 
for instance, be divided into two; one sub-catalog with 
clouds abiding in ambient tidal fields less than some ar- 
bitrary tide Tn m , and one with tide greater. The analy- 
sis of sub-catalogs into distribution functions can then be 
undertaken. It was found that clouds in low tidal field 
environments display significantly different characteristics 
than those in high tidal field environments. 

When represented in a log- log plot, the equivalent width 
distribution function (EWDF) is accurately fit by a power- 
law. The left-hand illustration of Fig. ^ shows the EWDF 
of Lya clouds, using all data irrespective of tidal field 
(heavy solid line). The break in the distribution, seen 
at log(W) ~ 1.5 (W ~ 32 mA, where W is the rest equiv- 
alent width of the Lya absorber in mA) is suggestive of 
two distinct populations of clouds represented in the data. 
The thin solid and dashed lines represent number weighted 
fits to the data at W < 32 mA, and W > 32 mA, respec- 
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tively. The values of the slopes are shown. These results 
strongly suggest that isolated clouds exist in a different 
kind of environment than un-isolated clouds. 
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Fig. 1. — The left-hand figure shows the log of the mean cu- 
mulative line density (solid jagged line) as a function of the log of 
the EW for low-z Lya clouds. It is evident that there is a broken 
power law, suggesting the existence of two populations of clouds. 
The spectral slopes, S, are shown for weighted linear fits to the 
mean EWDF for W > 32 mA (dashed line), and < 32 mA (thin 
solid line). The latter fit matches the void EWDF slope well (dot- 
ted line). The vertical dotted line is at 32 mA. The right-hand 
figure shows the trends of fitting parameters, based on the model, 
log(d\T fdz) = C + .Slog(W/63mA), with tidal field T Hm for void, 
and non-void catalogs. The upper panel of this figure shows the 
slopes, and the lower, the intercepts (log of line density at 63 mA). 
The horizontal dotted lines show the slope and intercept of the mean 
EWDF. Upper and lower sub-panels show non-void (tide as lower 
limit) and void (upper limit) EWDFs, respectively. The three ver- 
tical lines show the apparent range and center of the "transition 
zone" (see text). Note that characteristic slopes of void EWDFs 
(right-hand plot) are in agreement with the low EW slope of the 
mean EWDF (left-hand plot). 

Note that the EWDF of isolated clouds (Fig. □ left- 
hand plot, dotted line) closely matches the fit to the low 
EW part of the mean EWDF. This convincingly shows 
that isolated clouds constitute the "second" population 
(the low-EW side) seen in the mean EWDF. As it happens, 
there are virtually no low EW clouds in the un-isolated 
space (see e.g., Fig. 13 in Paper 1). I parenthetically note 
that the line density of void clouds at Nhi — 10 14 cm~ 2 
(i.e., about 225 mA) is on order 7 per unit redshift. This 
is about 1/5 of the mean value; certainl y not orders of 
magni tude less, as implied by simulations of lRiediger et all 
(199j)- 

Henceforth, I shall identify voids as regions that are iso- 
lated, where isolation is determined by the smallness of 
the scalar value of the summed tidal field (see §3 of Paper 
I)- 

The line density of void clouds at W > 15.8 mA 
is dN/dz ~ 500, and the spectral slope is steep S = 
dlog(d 2 Af/dzdW)/dlogW £ -1.5 ± 0.06 (Paper 1, Ta- 
ble 2). In contrast, the slope of the volume-weighted (i.e., 
mean) EWDF for W > 32 mA is S ~ -0.9 ± 0.1, while 
non-void clouds have quite flat EWDFs; S ~ -0.49 ±0.12. 
Void clouds appear to actually dominate the mean EWDF 
at low EW. 

The right-hand plot of Fig. ^shows the trends of fitting 
parameters with changes in limiting values of tide Tu m , for 
both void (lower sub-panels) and non-void catalogs (upper 



sub-panels). The tide is given in units of inverse Hubble 
time squared (see Paper 1) so that values are of order 
unity. There is a relatively narrow transition zone sepa- 
rating clearly distinct (linearly varying) trends in fitting 
parameters for void (low T) , and non-void (high T) envi- 
ronments (note the vertical dotted lines). The transition 
zone provides a phenomenological basis for establishing the 
tidal field which separates void from non-void clouds. In 
this way, it is shown that the space con taining void-type 
clouds fills - 86± 5 n % of the universe fsee lManning I (2003) 
for more details). This i s in good agreement w ith the 90% 
filling factor implied bv lCen fc Ostrikerl (|1999f> . 

1 summarize these introductory comments by noting 
that (1) N-body/hydro simulations appear to predict that 
the great majority of low-z Lya clouds are in the diffuse fil- 
amentary structure around galaxies, and not in voids. (2) 
The environmentally-resolved observations of Lya clouds 
in Paper 1 show that most observed clouds lie in the un- 
shocked void region, even without corrections for the much 
smaller redshift coverage of spectra at high sensitivity 2 . 
(3) The lack of the ability to resolve cloud halos in N- 
body/hydro simulations may be the principal cause for 
the inability to predict the existence of this population of 
void clouds. 

The foregoing tentative conclusions already throw sus- 
picion on some of the results of the simulations. Is it that 
simulations fail to predict void clouds entirely? Two ques- 
tions need to be answered; first, can diffuse clouds, ex- 
panding with the Hubble flow, explain the observed void 
clouds? Secondly, if not, can sub-galactic halos provide 
the gravitational potential needed to retain enough gas to 
produce detectable absorption systems? This is, in short, 
the purpose of this paper. 

The paper is organized as follows: In £J2I I investigate 
whether diffuse sheets of gas are capable of explaining the 
observed void EWDF, as envisioned in simulations. I find 
that diffuse sheets expanding with the Hubble flow are not 
consistent with observations. The remainder of the paper 
assumes the clouds are associated with sub-galactic halos, 
and that baryons are subject to pressure gradients after 
the reionization of the universe, and tend to expand away 
from their halos. In <J|ll present a general analytical treat- 
ment of spherical clouds to show the relationship between 
the baryon density profile and the slope of the column den- 
sity spectrum. In 2]I introduce the hydro/gravity simula- 
tion that will chart the evaporation of baryons from dark 
matter halos, in response to the re-ionization of the uni- 
verse. These halos are initially viewed as fully formed and 
not accreting. In I present the results of these simula- 
tions, based on a range of halo circular velocities and two 
halo types. Section presents a series of strategies for fit- 
ting model column density spectra to the observations of 
Paper 1. In fQl allow that clouds may accrete sub-halos, 
and so grow continuously from z = 6.5 to the present; the 
products at z = are presented, and the results are com- 
pared to those with fixed halos. In JEII reconcile Doppler 
parameters of model absorbers with those of the observed 
void clouds. In <J5]l summarize my findings. 

2. CAN DIFFUSE VOID CLOUDS EXPLAIN OBSERVATIONS? 

2 "High sensitivity" means sensitivity to very low EWs. 
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To determine the nature of void clouds under the cur- 
rent paradigm, one needs to know the fraction of neutral 
hydrogen atoms. Assume that clouds are optically thin, 
and immersed in a radiation field w ith an averag e flux 
J v = Mvlvo)~" s , where a s ~ 1.8 l|Zheng et all Il997t 
Shu ll et all . |l999) is the spectral index shortward of the 
Lyman limit vq. For the local universe, I assume, 

J = 1.3 x 1(T 23 erg cm" 2 s" 1 Hz" 1 sr _1 , (1) 

( Shull 'eFaTI |l999). For ionization equilibrium, 



nni = n H n e - 



(2) 



ijOsterbrockl 11989ft . where is the case-A recombina- 
tion rate coefficient appropriate for diffuse clouds. The 
photoionization rate is given by, 



T HI = 47T 



*)^- 3.24X10- (£)s- 



(3) 

where J23 = Jo/10 23 , and <t„ ~ ^(vq/v) represents 
the hydrogen cross-section for absorption of a photon of 
frequency v < i/q, where era — 6.4 x 10~ 18 cm -2 . 
The expected neutral fraction is, from Eq. [21 



In, ^ M,7 / 



a H {T) 



(A) 



J23J \2.51 x 10" 



n H . (4) 



Note that in a highly ionized plasma, n e ~ n.jj(l + 2%), 
where \ — n He/riH — 0.079. Given that O m = 0.3, 
0(,/0 TO = 0.1, and Y„ = 0.24, the hydrogen density is 
tih = 5.745 x 10~ 7 h^ 5 Cl c i cm -2 , where fl c i is the total 
mass density of a cloud in units of the critical. Given 
fi TO = 0.3, the mean hydrogen density is 



n H = 1.725 x 10~ 7 h% 5 



cm 



(5) 



^■75 ^cl- 



For T = 20,000 K, f HI = 1.55 x 10~ 6 

I note that the time-scale for ionization equilibrium is 
small; 

1 1x10V 



eq 



Thi 



J23 



(G) 



Thus, ionization equilibrium is a safe assumption. 

Now assume that large relatively homogeneous clouds 
exist in voids, as generically predicted by simulations, and 
that they are producing the isolated absorbers discovered 
in Paper 1. Assume that fV/f2 m = 0.3 (where J7y is the 
total mass density in voids). From Eq. the neutral 
fraction is, 



Jhi = 9.38 x Wn H T 



-0.702 1 2 
n 75' 



(7) 



where the index of T is the slope of the recombination 
c oefficient llOsterbrockLll989|) . 

iDave fc T ripp ( 200J) derive a relation between the tem- 
perature and the density from simulated data, 



T = b000(Sp/p) K 



(8) 



The exponent, 0.6, is actually an increasing function of de- 
creasing redshift and increases with reionization redshift, 



but .6 appears consistent with z = values llHui &: Gnedinl 
Il997|) . According to Eq. 4 of IDave fc Trind lj2001|) 
Sp/p « 12(A^ ff/ /10 14 ) - 7 at z = 0. The over-density re- 
quired to produce a column density N HI — 10 12,5 cm -2 
is 0.379 n H , or 6.54 x 10~ 8 h^ 5 cm- 3 . With Eq. 01 de- 
rive a temperature T — 2794 K. Substituting this into 
Eq. |3 the neutral fraction of hydrogen is found to be 
Jhi — 8.50 x 10 -6 . Thus the length of the column needed 
to produce the absorber is, 



I = 



N H i N H i 



Shiuh 



1.85 Mpc. 



(9) 



If the filament is expanding with the Hubble flow, then 
the velocity width, ignoring thermal broadening, would be 
~ 140 kms . Because the background mass density in 
voids is expected to be about an order of magnitude less 
than in the filamentary medium surrounding galaxies, the 
frequency of fluctuations of amplitude Sp/~p = 0.379 will 
be very much smaller in voids than in non-void regions. 

For N H i = 10 13 cm- 2 , Sp/p ~ 1.07, f HI ~ 8.23 x 
10" 6 cm" 2 , T = 5207 K and I ~ 2.14 Mpc. The im- 
plication is clear; absorbers formed in this fashion would 
be seen only as very broad troughs. The Hubble flow 
across I would then be well beyond the values found in 
the Space Telescope I nfrared Spectrograph Echelle spectra 
ijDave fc Trippll2001l). In GHRS sp ectra with resolution is 
~ 19 kms" 1 ijPenton et alll200 0a). void clouds have Lya 
Doppler parameters, b ~ 30 kms -1 (Paper 1), probably 
large by a factor of ~ 2 du e to resolution pr oblems, and 
the presence of components l|Shull et al.ll2000|) . Thus void 
clouds have Doppler parameters only 1/5 to 1/10 as large 
as that which diffuse sheets would produce in voids. 

In conclusion, it appears clear that the detected void ab- 
sorbers must be either significantly denser than the mean 
density (i.e., p/~p w 10), or the velocity field of the baryons 
must be strongly convergent with respect to the Hubble 
flow, or both. Neither of these conditions is consistent 
with the model of void clouds as diffuse baryonic slabs 
expanding with the Hubble flow. Thus, observations ap- 
pear to be in conflict with the paradigm promoted by the 
N-body simulations as applied to Lya clouds. 

In the following section we begin our analysis of sub- 
galactic perturbations as a parent population for void 
clouds. 

3. AN ANALYTIC TREATMENT OF SPHERICAL CLOUDS 

If homogeneous slabs of gas expanding with the uni- 
verse cannot explain these clouds, then what might? I 
hypothesize that absorbers stem from sub-galactic pertur- 
bations that were massive enough to restrain the outflow 
of gas from the dispersive tendencies of the Hubble flow 
and the pressure gradients that result from the ionizing 
background. I start with an analytic treat ment of clou ds 
taken as baryons held within a "mini- halo" lReellll986^ . 

Clouds whose evaporation is restrained by a centrally 
condensed halo of dark matter will have a radial density 
gradient. If such clouds can be roughly approximated by 
power law density profiles, their ionization structure, and 
column density spectra can be calculated. 

3.1. The H I column density of spherical clouds 
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I initially assume that the baryons of a cloud have a 
spherically symmetric power-law distribution with index 

p(r)ocr-' 5 . (10) 

Eq. 01 shows that uhi oc n 2 H . The H I column density 
produced at projected radius r p from the cloud center is, 



N H i 



n HI (r)dl, 



(11) 



where I is constrained by the equation r 2 = I 2 + r 2 . From 
Eas. EH and ITU we find. 



2a 



(^) 



*-(l+2 X )/ (n H {r)fdl. (12) 

1 HJ JO 



3.2. The Column Density Spectrum of Spherical Clouds 

What column density spectrum would a cosmological 
distribution of such clouds produce? The cumulative col- 
umn density spectrum may be broken down by the chain 
rule, 

!(Nm) _ m^sii _ «w * (13) 



92 97V, 



<9z 9r cW 



HI 



where AT = n c i(> N^j). The relative number A/" of ab- 
sorption systems produced by a spherically symmetric sys- 
tem with a column density > Njjj at projected radiu s r p 
has a derivative, dM/dr p oc 27rr p . By Eqs. l2l and ITOl 



uhi oc r 



-2# 



Then by Eq. ITTI 



Nmcxr-W-V. 



(14) 



Carrying out the derivatives in Eq. 1131 and simplifying, 

f{Nm) -dN H l{^) <xNHI ' (15) 

which is the differential distribution. 

The cumulative column density distribution has a slope 
larger by 1; 



S = 



1-20' 



The inversion of this yields, 



5 2 



(16) 



(17) 



These two equations provide an important diagnostic of 
the observed column density spectrum (CDS) if we are to 
interpret it as produced by discrete, spherically distributed 
clouds, because it determines the cloud density profile that 
is consistent with a given CDS slope. Note that the steeper 
the index of the density profile the flatter the spectral 
slope S. For instance, the observed slope of the void cloud 
EWDF, S ~ —1.5 can be explained if ■& ~ 7/6; i.e., p oc 
r -117 , while the slope for non void clouds (S ~ —0.5) is 
explained if p oc r~ 2 5 . 



4. NUMERICAL SIMULATION OF CLOUDS 

For baryons associated with dark halos, it is not ex- 
pected that baryons woul d follow their halo mass distri- 
bution at late times (e.g., iBarkana &: Loebl H"999(l - one 
might expect a flatter distribution owing to the effects of 
evaporation. As noted above, the steep EWDF of void 
clouds (Fig. rij imply a relatively flat density profile; this 
may be a measure of the effects of reionization on baryons 
held by small dark matter (DM) halos. These consider- 
ations suggest that the clouds producing the absorption 
in voids have evolved into a rather flat density profile. A 
more accurate view of the probable profiles of such clouds 
may be achieved by a cosmological simulation of the re- 
sponse of DM-held clouds to reionization. 

I propose to follow in detail the likely evolution of 
baryons in clouds grav itationally held by dark halos before 
reionization (z ~ 6.5. iDiorgovski et all l200i|) . As noted 
above and in Paper 1, 3-D simulations suffer from space- 
and mass-resolution problems when small-scale structures 
are of interest, making it difficult to separate actual phys- 
ical results from numerical artifacts. When small-scale 
structure needs to be studied, a 1-D code may provide a so- 
lution to the problem. I acquired a copy of a 1-dimensional 
Lagrangian hydro/gravity code o flThoul fc Weinberd l)1995f) 
courtesy of one of its authors, (A. T.), and adapted 
it to the needs presented by the problem posed above. 
Among the adaptations made was the transformation from 
Einstein-de Sitter universe to a flat FRW (A = 0.7) with 
h = 0.75, and the writing of initializing routines for vari- 
ous halo models. These are consistent with the parameters 
used in Paper 1. A cosmic baryon ratio £lb/Q m = 0.1 is 
assumed. 

4.1. Nature of 1-D Code 

The code sets the value of the gravitational constant G 
to unity. In addition, it sets a mass scale, M = 1.0 x 
10 11 M Q , and a length scale L = 1.0 X 10 12 AU, or 22.592 
kpc, mandating a time scale t — 5.042 x 10 15 s, or ~ 160 
Myr. 

4.1.1. Governing physics 

The equations, and related details of the computational 
scheme that forms the basis for tracking the evolutio n 
of the clouds are detailed in iThoul fc Weinberd l|1995|) . 
The basic components are recapitulated as follows. Ion- 
ization equilibrium is assumed because its time-scale is 
smaller than the time-scales of interest. The energetic 
gain r from photoionization of H I, He I and He II by 
UV photons is calculated on the basis of self-consistent 
concentrations of these components. We assume that the 
spectral index above the Lyman limit foll ows the rela- 
tion, J v = Jvi{i> I Vft)~ a , where a = 1.8, (e.g.. IZheng et all 
119971 IShull et all 119991) . The redshift var iation of J (z) 
is adapted from IThoul fc Weinberd l)1995[) . and is shown 
graphically in Fig. [2J The radiative losses, A, are found 
by summing the cooling rates from collisional excitation 
and ionization, recombination, and br emsstrahlung, given 
the io nic abundances and the density ijThoul fc Weinberg 
1995). In the La grangian code, bin boundaries for baryons 
respond to the effects of pressure differentials and gravi- 
tational potential gradients. The change in the internal 
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energy per unit mass u with time is then, 



du 
~dt 



P 

Pi dt 



+ -(r-A), 

Pa 



(18) 



where the first term on the right hand side is the adia- 
batic loss; p is the pressure, and the g-subscript refers to 
corresponding values for gas-bins. Shocks occur when mo- 
tions are convergent. The a rtificial viscosity technique of 
Ric htmver fc Mortonl l)1967|) is used, adjusted to spread 
shocks over 4 or 5 shells, helping to ensure energy conser- 
vation l)Thoul fc Weinberg . 1 19951) . The pressure is given 
by the equation of state, 



Pgl 



(19) 



The baryonic component responds to forces as follows, 



dVg 

dt 



= -Anr„ 



dp M(r g ) 



dm„ 



(20) 



where M(r g ) includes dark and baryonic mass, dp is the 
pressure differential over a Lagrangian bin, and dm g is the 
gas mass contained therein. 
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Fig. 2. — Th e prescribed average inten sity at the Lyman limit, 
as adapted from Thoul & Weinberg 1 1995). In the interval < z < 
2.0, Jq(z) oc Jo(1 + z) 3 . The final results of simulations are fairly 
insensitive to the details of reionization. 

The time-step for the simulation is determined by the 
smallest of a collection of time scales, including the dynam- 
ical timescale, the Courant, and the cooling time scales, 
and a time-scale to ensure that fluid shells do not cross. 
In addition there is a maximum length time-step of 10~ 3 t. 

The minimum time-scale is frequently determined by 
the dynamical time-scale in the first few shells. When the 
timescale for cooling is much less than the dynamical time 
scale, shells are essentially in free-fall. This causes the cen- 
tral density to increase dramatically so that cooling time 
steps become extremely small. In response to this prob- 
lem, when t coo i <C Td yn the program "freezes" the shells 
in question, and refers the mass to a minimum radius, 
here, r m i n = r c — essentially taking the mass in the shell 
out of gas form (we might speculatively suggest that gas 



fragments and forms stars in the collapse zone). In prac- 
tice, this does not occur for syste ms of circular v e locity 
less than v c <J 21 or 24 kms" 1 , for iNavarro et all ([l996, 
1997) (hereafter, NFW) and isothermal halos, respectively. 
Since this is an investigation of clouds and not galaxies, I 
ignore systems of larger circular velocity. I check that the 
point at which shells are frozen is physically realistic by 
monitoring the ratio of the diameter of the first Lagrangian 
bin to the Jeans length, given by 



GpfiUlH ' 



(21) 



where 7 is taken to be 5/3, and p is the mass per particle 
in units of the hydrogen atom, that approaches the value 
p = 0.588 at the low densities of model clouds at low to 
moderate redshifts. This indicator is found to approach 
unity shortly before the shells are frozen, indicating that 
the program is treating the problem correctly. 

It is possible that starting with a less-relaxed halo, 
warmed by secondary infall at z < 6.5, halos of circular ve- 
locity high enough to collapse in this way under the current 
arrangement, would not collapse. This implies that clouds 
of somewhat larger circular velocity could exist without 
forming stars (see 

4.1.2. General Treatment of Dark Halos 

Because the iThoul fc Weinberg! |l995) code was origi- 
nally written for collapsing protogalactic clouds, the pro- 
gram's treatment of DM is inappropriate for the present 
simulation. Instead, a dark halo overdensity is established 
and not allowed to change (these halos are to be specified 
below). It would be more proper to project a dark halo as 
it would be at a high redshift, growing somewhat by sec- 
ondary infall, so that the end product is consistent with 
the specified circular velocity. However, of this there are 
two views; one state s that the halo of a s mall cloud is in 
place by z — 6.5 f see IWechsler et aT1l2002|) 3 . so that most 
of the response to absorbed radiant energy is at redshifts in 
which the total cloud mass is fairly stable. In this scenario, 
therefore, assuming a fixed halo circular velocity can only 
modestly minimize the final extent of cloud expansion due 
to its greater average potential. However, the other view 
is based on observations suggesting the growth of halos, 
even at late times. Therefore, after applying the analy- 
ses to fixed halos, we then allow halos to grow steadily 
by accretion of sub-halos, leading to a growth of circular 
velocity with time. These "grown" halos are subjected to 
the same analysis as the fixed halos, providing a plausible 
maximal range for reasonable models of halo growth. 

The spherically symmetric conditions of the simulation 
presume that the halos are isolated from the influences of 
others. This approximation should be accurate for void 
clouds since the environment around them is expanding 
with the Hubble flow, and evaporation ejecta are not likely 
to interact with that of other isolated systems. However, 
it would not likely apply to non-void space because in the 
convergent flows of over-dense regions, cloud outflows may 
collide causing shocks, that violate the isolation and spher- 
ical symmetry assumed in this simulation. 

3 The virialized mass will grow in time, but that is only apparent 
growth since it results from the decline in p cr it with time. 
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4.1.3. General Initial Conditions 

The initial conditions of model clouds are designed to 
make them dynamically stable. Because the central densi- 
ties of clouds are supposed to approximate the density of 
the universe at the redshift of their formation (the canon- 
ical explanation for the high concentration parameters of 
small systems in relation to those of large systems found 
in simulated halos), I assume that the halos I model are 
already in place at z = 6.5. There may be some confu- 
sion on this issue because it is seen that the virial mass of 
even small ha l os co ntinues to grow down to z = (e.g. 
iBullock et all l2001j) . However, the virial radius is ap- 
proximately equal to the radius within which the average 
density is ~ 200 p cr it- Therefore, the decline in the back- 
ground density with the expansion of the universe causes 
an expansion of the virial radius, resulting in an increase in 
the virial mass with time. Comparisons with code kindly 
lent from J. Bullock ("cvir.f") show that, for clouds with 
v c < 25 kms" 1 , there is very little increase in mass, or 
decrease in the concentration parameter over and above 
that which naturally results from the expansion of the uni- 
verse. Therefore, from the standpoint of current thinking, 
assuming that the halos were in place at z = 6.5 should 
be non-controversial. However, from the standpoint of the 
isothermal model, this may not be true. Accordingly, we 
do investigate halos which grow steadily in §J| 

In deciding the initial conditions of the gas associated 
with a "fixed" dark halo (i.e., fixed circular velocity v c and 
overdensity profile), it is important to note that at high 
redshift the baryonic cloud is immersed in a critical density 
background with density comparable to that of the cloud. 
For, at z=4.84 the background density is 200 times the 
current critical density, thereby making the virial radius 
significantly smaller than at the present. Because baryons 
are dissipational, there is no way to separate the baryons 
that are part of the overdensity from those that are part 
of the background. In fact, what allows for the existence 
of a prolonged and stable overdensity is a distortion in the 
Hubble flow that is concentric about a point we would call 
the center of the cloud. Since these halos are immersed in 
a high-redshift universe, their overdensity exists by virtue 
of this deviation from the Hubble flow. 

At the initialization stage, an outward motion is given 
to the gas that is less than, but gradually approaches the 
Hubble flow velocity as one considers radii farther from 
the cloud center in a manner such that the overdensity 
is stable in time. This is actually easier than it sounds, 
for one need only give an outward velocity such that the 
radial flux of gas is what it would be if the universe were 
flat and of critical density. Thus, at a given redshift z, the 
outward motion of gaseous matter v g needed to maintain 
a density p g (r) at a distance r, for instance, should equal, 

v g (i) = Hr(i)n b p crlt /p g (i), (22) 

where the variables, H, Qf,, and p cr it are evaluated at the 
redshift in question. This creates an overdensity that is 
stable within the Hubble flow (the peculiar velocity sup- 
ports the overdensity) . The dark matter is given the same 
peculiar velocity field. The velocity field that results can 
be expressed as proper radial motion (including the Hub- 
ble flow) as well as peculiar velocity, with the Hubble flow 
subtracted. 



Because of the exigencies of a Lagrangian code, and 
the apparent lack of evidence for a truncation of halos 
at the virial radius, it was decided not to impose a radial 
truncation of the overdensity of halo models. The halos 
were extended to a distance at which the average inter- 
nal overdensity in units of the current critical density was 
p = 2 p cr it , equivalent to a distance of ~ 2 Mpc from an C* 
galaxy with an isothermal (hereafter referred to as ISOT) 
model halo. Later, when estimates of the average density 
contributed by clouds are required, it will be possible to 
impose a truncation radius, such as that featured in Pa- 
per 1. It should be borne in mind, however, that there are 
no features in the simulation products that would imply 
that the mass is truncated at any characteristic radius, nor 
does the preponderance of observational evidence (see be- 
low) require it; only N-body simulations suggest something 
like that might be required. Since the extent of gravity is 
infinite, halos have the tendency to draw matter toward 
them incrementally at great distances. 

Initially, baryons are taken to "follow" the DM, so there 
is an inward peculiar gravity acting on the baryons. To 
restrain baryon collapse, I initially constrain the cloud to 
have a temperature such that the pressure exactly balances 
the peculiar gravity in each bin. This will not guarantee 
that clouds are static because cooling rates vary within 
the cloud, however, it is dynamically stable. For instance, 
if the simulation is started at z = 10, instead of 6.5, the 
baryons will "slosh" about within their halo, but when 
ionization begins, the cloud expands, producing a result 
remarkably similar to those with simulations begun at z = 
6.5. 

It turns out that the details of the temperature struc- 
ture is relatively unimportant because, for clouds with 
virial temperature T V i r < few x 10 4 K, the cooling time- 
scale is large, suggesting that clouds are fairly stable before 
reionization. Further, the energy deposited upon reioniza- 
tion dwarfs the initial energy content, so that the effect of 
the initial temperature on cloud dynamics is insignificant. 
Other tests showed that the final state of the cloud is also 
quite insensitive to variations in the specific prescription 
for the onset of the ionizing flux (see Fig. |5J). Thus, final 
results are sensitive only to halo type and circular velocity. 

4.1.4. Initializing routines 

In the initialization stage of the program, the halo model 
and its DM circular velocity are specified, and the total 
mass of the cloud in DM and baryons derived therefrom 
are calculated according to the cloud halo model, extend- 
ing to a radius at which the internal over-density is 2 times 
the critical. This ensures a large enough scale to capture 
the relevant dynamics of expansion, and guarantee that all 
observable column densities occur within the cloud end- 
products. The mass is then divided equally into the num- 
ber of Lagrangian bins to be used (generally, the number 
of gas bins are set to n g — 200, and the dark matter di- 
vided into nd = 500 bins). Then the initial radius at which 
these masses are found is determined by straightforward 
application of model relations (see §5.1). 

This study considers halos which are isothermal as well 
as NFW halos. I do not truncate these halos at the virial 
radius because what is of interest is the mass distribu- 
tion of the halos, rather their kinetic properties. Nu- 
merous studies now support the claim that there is no 
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sign of a drop in the index of the density profile be- 
yond the virial radius, or for that ma tter, the NFW ra- 
dius of maximum circular ve locity (s.^.- I Zari^ sky Whi^d 
199 4 IZaritskv et all Il997t IZabludoff k. Mulchaevl 12000: 
McKay et all 120021) . In their study of weak lensing, 
McKav et all l)2002|) show that the line-of-sight (LOS) ve- 
locity dispersion ay of satellite galaxies around 618 iso- 
lated host galaxies from the Sloan Digital Sky Survey is 
consistent with ay being constant over a range 100 <J r <J 
500 h _1 kpc, and that the number density of satellites 
n oc r~ 21 over this range. This is very suggestive of a 
mass distribution that is not truncated; one not inconsis- 
tent with an isothermal mass distribution. To make the 
distinction clear, when I refer to a halo, I am referring to 
the mass distribution, rather than the mass in the viri- 
alized portion thereof. It is expected that the mass of 
the early halo will generate peculiar velocities in the sur- 
rounding universe, contributing to an enhancement of the 
density outside the virialized region of the halo. 

The NFW Halo: The NFW halo is treated in consider- 
able detail in Paper 1. Equations 34 to 40 of Paper 1 de- 
scribe the derivation of the scaling relations of the model, 
which are necessary for calculating halo parameters as a 
function of halo circular velocity. A formula for the mass 
as a function of distance is the basis for determining the 
bin distributions and calculating density within the bins, 

, s 4tt , | log e (l + r/r a )- jgf- \ 

M NFwW = T 200p c ^3 00 , 

(23) 

where r s is the scale length, r2oo is the assumed virial 
radius enclosing a density of 200p C Ht, and c is the concen- 
tration parameter defined by, c = r2oo/ r s, where r s is the 
scaling radius, and r2oo is the density at which the aver- 
age internal density is 200 times the critical density, or, 
loosely, the virial radius. The concentration parameter is 
a function of the halo circular velocity, and was derived 
from relations taken from NFW (see §5.1.2 of Paper 1), 

c~ 132.3 «-° a f 3 > (24) 

where v max is the maximum halo velocity, which occurs 
at a distance, r rnax ~ 2.16r s . Be aware that there is sig- 
nificant dispersion in these relations, as they are average 
values derived from the results of computer simulations. 
The NFW halo has a cuspy core, with density approach- 
ing p oc r" 1 as r goes to zero. 

The Isothermal Halo: In Paper 1, the isothermal halo 
(which I will henceforth abbreviate as ISOT) was used to 
model the total mass of galaxies, and so required no inter- 
nal structure to be stated. This study, however, requires 
an explicit core. I assume that the central non-singular 
co re of the ISOT model has a density equal to that found 
by iFirmani et all l)2001(l in their investigation of the cen- 
tral densities of dwarf galaxies. I apply an exponential 
term to the density law, 

p(r) = ^(l-e-C^') 1 (25) 

so that that there is a flat core of density p c ~ 0.02 M Q pc _ 
that constrains the core size ro, and also allows the mass 



to approach a value M(r) — Kf for r ^> ro. I find that 

r ° = °- 3035 x (idb*) kpc - (26) 

The simulation begins with the baryons "following" the 
DM distribution, and shortly thereafter they are subjected 
to an isotropic metagalactic flux whose amplitude at the 
Lyman limit Jq(z) is shown in Fig. [21 
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Fig. 3. — The evolution of an isothermal cloud (left panel), and 
an NFW cloud (right panel), both with v c = 21.1 kms -1 . The log 
of the baryon number density is shown as a function of the log of the 
cloud-centric radius at an arbitrary series of times (see text). The 
simulation starts at z = 6.5 (topmost lines) and proceeds to z = 0. 
The NFW halo (right-hand plot), being more strongly concentrated, 
has a smaller first gas bin, and higher density there, but a lower total 
mass. Note the outflowing wave that is stronger in the NFW halo 
than the ISOT. 
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Fig. 4. — The proper velocity (upper panels) and the peculiar 
velocity (lower panels) for ISOT halos (left panels) and NFW halos 
(right panels) at arbitrary points from z = 6.5 (left-most line) to 
z = (right). As in Fig. |3] velocities are shown as a function 
of the log of the system-centric radius. The cloud halo velocity 
shown in this experiment is Vc = 21.1 kms -1 . This NFW halo has 
acoustic oscillations on time-scales of order 600 Myr. The ISOT 
model produces relatively minor shocks and oscillations. 
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4.2. Behavior of Baryon Clouds after Reionization 

The program advances in steps and needs about 15,000 
time steps to get to z — 0. The program prints out data 
every 20 steps, leaving about separate 750 files. Figures 
and0]show the reaction of a large cloud (v c ~ 21.1 kms -1 ) 
to reionization. The time steps at which the distributions 
were sampled are arbitrary; chosen to give a sense of the 
dynamics and to avoid confusion. The chosen files are 
roughly geometrically arranged {e.g., using files #1, 2, 5, 
10,..) until (for ISOT) reaching file number one hundred, 
then continued to sample at each 100 steps to the last step 
near 750. For NFW halos the doubling continues to z = 
to avoid the confusion of overlapping plots. A sense of the 
rcdshift can be had by recognizing that the density at the 
extreme edge of the cloud is flmPcrit/mH oc (1 + z) 3 . 

The principal distinction between the ISOT and NFW 
density profile evolution is that ISOT has a softer out- 
flowing wave, and a significantly larger central density at 
the conclusion of the run. The NFW halo, which is very 
concentrated, sends out a strong wave of expansion that 
introduces a shock that reverberates back into the inner 
cloud. The lower final central baryon densities of the NFW 
clouds, roughly a tenth of the ISOT clouds when com- 
paring clouds of equal circular velocity, reflects its lower 
mass. The evolution of velocity profiles for halo models 
with v c — 21 kms -1 is shown in Fig. It shows that 
the outflowing wave from the NFW halo has a signifi- 
cantly stronger amplitude - by almost a factor of two - 
accompanied by a strong oscillation of the inner cloud. 
The magnitude of the wave declines significantly in time 
for each model. The final peculiar velocities (lower pan- 
els) at large cloud-centric distances are only slightly nega- 
tive for the NFW halo, while the more massive ISOT halo 
has produced a sizable negative peculiar velocity at late 
times just beyond the outflowing wave. Note, however, 
that the higher density gas internal to this still has a pos- 
itive peculiar velocity. This discontinuity, here located at 
r ~ 260 kpc, is far beyond the canonical virial radius for 
this cloud (r V i r ~ 25 kpc). 

4.3. Comparison of Program Output with Other Work 

Data from simulations of clouds can be used to probe 
the trend of temperature with density as a function of 
rcdshift. In Fig. [SJis plotted the logarithm of the density 
(in units of f2t> Pcrit) versus the log of the temperature 
at z=0 (skeletal), z = 1 (open circle), and 2 (open pen- 
tagon). The data at z=0 is taken from a maximal ISOT 
halo with v c — 23.7 kms -1 (5-pointed stars) and with a 
smaller halo, v c = 13.3 kms -1 (crosses). Note that the 
locus of points from the smaller cloud is marginally cooler 
than that from the larger cloud because the smaller cloud 
may expand more freely, increasing the adiabatic cooling. 
The heating produced by the outflowing shock (or cooling 
following its passage) can be readily seen in the z = 1 and 
2 plots (note the discontinuity) , though by z — the wave 
has passed beyond the outermost point. 

Also in Fig. [5] these pro gram outputs are co mpared 
with a result from Fig. 36 of lHui fc Cxnedinl lHool . based 
on the output of a 3-D hydro code (dashed line, label 
HG, z = 2), which utilizes rapid reionization at z = 7 
and a flat A model. A lso shown is the rela tion implied 
by Eq. 03 as found in iDave fc Trippl l)2001|) (solid line, 



labeled DT), which represents the trends in the filamen- 
tary structure surrounding galaxies (i.e., n o n- void space) 
at z ~ 0. Note that the iHui fc GnedirJ (|1997j) trends 
(z = 2) appear to approach the trend for the results from 
the 1-D simulation insid e the outflowing wave. However, 
the|D ave fc Triod J2001 model appears to predict signif- 
icantly larger temperatures at densities of a few or more. 
This is attributed to shock heating and bulk motions in 
this denser environment, said to be more pronounced at 
low redshift lDave fc Trippl (|2001^ . 

Given the differing physical conditions that are probed 
by the 3-D (non-void) and the 1-D (void) simulations, the 
differences in predicted temperatures are not surprising. 




Log(l + 5) 

Fig. 5. — The logarithm of the temperature (in Kelvins) of 
cloud parcels as a function of the log of the over-density (in units of 
f^b Pcrit) for various redshifts. Symbols represent the p r esent work, 
and labeled lines rep resent the work of|Hui & Gnedin (1997), and 
Dav e fc Trippl 120011) . The skeletal symbols show the variation of 
trends of temperature at z = in large clouds (5-armed) vs. small 
clouds (4-armed), where smaller clouds are cooler as a result of 
the greater freedom a small cloud has to expand, and adiabatically 
cool. Open circles, and pentagons show trends for z = 1, and 2, 
respectively, for the large (v c = 23.7 kms -1 ) cloud. See 14.31 for 
further discussion. 

5. RESULTS AT REDSHIFT ZERO 

The program is run for a range of cloud halo velocities 
v c with bin-centers at, 

Vc (n) = 5.31 x 10 - 05 ™ kms" 1 , (27) 

for n=0 to 13. These bins cover circular halo velocities 
of 5 to 25 kms -1 . Halos larger than this experience what 
appears to be the equivalent of star formation, i.e., free-fall 
collapse in the central regions of the cloud. The cadence 
in Eq. 1271 has 20 logarithmic steps in a decade of v c . 

Another form for designating the velocity bins is the 
following, 
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One should notice that 144 kms -1 is the velocity of the 
dark halo of an C* galaxy that, when baryons are added 
to comprise 10% of the total, has a circular velocity of 161 
kms -1 (see §5.1.1 of Paper 1). 
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Fig. 6. — Final density profiles (z = 0) for the range of halo 
velocities for the ISOT model. The dashed line represents the slope 
of the density profile = 1.125) necessary to pr oduc e the observed 
slope S RJ — 1.6 of the CDS for void clouds, (see Hi}. II . Clouds with 
circular velocity v c > 10 kms -1 have slopes that approximate this 
va lue inside ~ 150 kpc. Profiles with dotted lines are undetectable 
in IPenton"*c t al. (2000aj). The enlarged dots and short-dashed lines 
represent the radii where particular clouds give the stated column 
density (see log of column density labels). The inset shows the 
detail of the cloud outer edge showing that the self-gravity of the 
larger clouds has restrained the outflow wave for halos with velocity 
v c > 10 kms -1 . 

5.1. Model Halos at z = 

The products of simulations for z = are of special in- 
terest, for the density profiles, temperatures, and neutral 
fractions of hydrogen, will provide what we need to pro- 
duce a model column density spectrum (CDS). The results 
strongly depend on the halo model. 

5.1.1. ISOT density profiles 

The z = results of an ISOT simulation are seen in 
Fig. El where the smallest cloud v c = 5.31 kms -1 is the 
bottom, dotted line. The smallest cloud to produce a de- 
tectable absorption line of Njji > 10 12 ' 4 cm -2 is the bot- 
tom solid line, wherew c = 10.6 kms -1 . The long-dashed 
line in the figure has a slope d that is associated with a 
CDS slope S — — 1.6 through Eq. ^1 the observed slope of 
the void CDS. The labeled short-dashed lines show the ra- 
dius at which selected column densities are detected, and 
the baryon density at the closest approach of the LOS. 
This will be discussed in detail in ^5.31 Note that large 
clouds may produce low-column density absorbers at dis- 
tances as large as 200 or 300 kpc, though the spectral slope 
they would produce would be significantly shallower than 
that observed in the void CDS of Fig. ^ 



The inset of Fig. shows an enlargement of the zone 
of the outflowing wave. For smaller halo velocities, the 
position of the wave at z = scales directly with the core- 
radius, indicating the outflow is essentially unrestrained. 
However for v c ;> 10.6 kms -1 (8th from the top), the 
position of the wave is at progressively smaller radii, at- 
tributable to the cumulative effects of their massive halos. 
These effects are also seen in the peculiar velocity plot of 
the large ISOT halo in Fig. H 

That neutral hydrogen mass of these halos can be 
straightforwardly calculated from the output of the code, 
and while not directly relevant to this study, it is relevant 
to 21 cm studies of high velocity clouds (HVCs), and is an 
opportunity to compare Lya and HVCs. These data are 
presented in <J7| 

5.1.2. NFW density profiles 

The NFW simulations to z = are presented in Fig. 
[7| Gas density profiles are seen to be flatter than § — 
1.125 outside ~ 30 kpc. There is an apparent total blow- 
out phenomenon for clouds with v max < 10 kms -1 , and 
a blowout of the outer cloud for clouds smaller than ~ 18 
kms -1 . The smaller physical extent of the clouds has 
resulted in the expanding wave having traveled beyond 
the outermost Lagrangian bin except for the most massive 
cloud with v c = 21.1 kms -1 . The cloud with v c = 23.7 
kms -1 , however, collapsed from the inside out, and would 
presumably form stars, and was therefore discarded as a 
model cloud. As we shall see, NFW halos do not produce 
observable H I columns in these simulations. 
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Fig. 7. — NFW density profiles evolved to z = for v c = 5.3 to 
23.7 kms -1 . Note that low v c clouds have expanded with such force 
that they are actually underdensities at z = 0. The dashed line has 
slope -& = 1.125, that would produce a column density spectrum 
slope S = —1.6, approximately that of the void CDS. 
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Here I discuss the method of extracting model absorp- 
tion systems, and looking at the thermal and velocity 
structure of the absorbers. 

5.2.1. Calculating cloud column density 

Let us consider a line of sight through a particular ISOT 
cloud. Given an arbitrary targeted column density, we 
seek to find two consecutive bin radii that serve as impact 
parameters r p (i) and r p (i + 1), at which column densities 
bracket the target column. One may then interpolate the 
solution. When the line of sight (LOS) velocity profiles of 
the model absorbers are analyzed, it is found that absorp- 
tion systems typically have a strongly peaked central core, 
with broad wings that may add substantially to the calcu- 
lated column density. Thi s "background" is attrib utable 
to a Gunn-Peterson effect l)Gunn & Petersorlll965|) . Sim- 
ulated clouds are hypothesized to be immersed in a gas of 
density Qbp cr i t , and under the UV flux utilized in these 
simulations (see Fig. [21 and Eq. ^) , the neutral fraction 
of hydrogen gas of this density at z — is found to be 
f HI ~ 1.5 X 1(T 6 cm- 3 . 

The simulated data is subjected to the same treatment 
as the observations. In the reduction of the observational 
data, a continuum is placed with respect to which the fluc- 
tuations of the spectrum, in zones not influenced by ab- 
sorption or emission lines, are evenly balanced, above and 
below the continuum. Therefore, this background Gunn- 
Peterson effect is subtracted from the model absorber. Us- 
ing the neutral fraction of this background gas, and the 
z = background hydrogen density (1.745 x 10~ 7 cm -3 ), 
yields a neutral hydrogen density in a cloud of total density 
QmPcrit of uhi = 2.61 x 10~ 13 cm~ 3 . The neutral column 
in a megaparsec is thus Njji = 8.04 x 10 11 cm -2 . In the 
Hubble flow, if we take 1 Mpc to be 75 kms -1 , there is a 
column "dispersion" of dNni/dv — 1.07 x 10 10 cm -2 per 
kms -1 , or 10 12 cm -2 per 93.3 krns" 1 . This provides 
a minimal approximation to the displacement of the con- 
tinuum due to background neutral hydrogen. This back- 
ground column density dispersion is scaled to the velocity 
of the extracted profile, and subtracted from the total col- 
umn density to get an "observed" value. 

The program finds the projected radius r p at which a 
target column density occurs, and uses that to calculate 
the surface area of the cloud that will produce a column 
greater or equal to the target column. The utility of this 
parameter for deriving the column density spectrum will 
be demonstrated in 



5.2.2. Cloud temperature versus b-value 

The temperature of a cloud is an important factor in 
predicting the Doppler parameters of model clouds. As 
shown in Fig. \5\ the temperature of clouds vary system- 
atically with density. I define the temperature of an ab- 
sorption line as the temperature of its core. To do this I 
calculate the column density- weighted temperature within 
the range at which the differential column in the successive 
bins drops to half the maximum value. 

Cloud outflows at z = are generally quite small (of 
order ~ 2 to 7 kms -1 ) compared to the observed b- values 
of clouds. Recall that the most likely Do ppler parame- 
ter for void clouds with the GHRS data of lPenton et alJ 
(20003) (see Paper 1, Fig. 6) is b Lya - 30 - 35 kms" 1 . 
If purely thermal, this would correspond to a temperature 



of T ~ 70, 000 K. I retrieve from the model lines a col- 
umn density-weighted temperature for a range of columns 
Nhi, and a range of cloud velocities, all summed within 
the line FWHM. This will provide a theoretical lower limit 
to observed cloud b- values as a function of EW or column 
density. In fact, a fairly tight relationship between column 
density and temperature results from model clouds. These 
temperatures, for log-column densities of 12.5, 13.0, 13.5, 
and 14.0 are, 7800 K, 11600 K, 16900 K, and 27000 K, re- 
spectively. The range of temperatures at log-column 12.5 
is fairly large (~ 300 K) because of a wide range of halo 
velocities that may produce the line, but it is smaller for 
large column densities since only a small range in cloud v c 
may produce them. These data, transformed into a ther- 
mal b- value 4 , are plotted in Fig. |S]as open pentagons, and 
provide a theoretical lower limit to cloud b- values.. 

5.3. Transforming EW into Column Density 

In order for observed absorption systems (EW in units 
of mA) to be compared with model results (column den- 
sity in units of cm -2 ), one needs to account for the self- 
blanketing effects of H I clouds. In this section I discuss 
transformations of the EWDF to a CDS, for this is re- 
quired before model clouds may be compared to the ob- 
served EWDF. 

If the cloud is quiescent and of uniform temperature, 
then the thermal broadening of the line will character- 
ize the Doppler parameter b. A low-temperature cloud 
will deplete the continuum photons with wavelengths A = 
1215.67 A in the cloud rest frame because of the large 
fraction of neutral atoms that can efficiently absorb at the 
Lya resonance, and so produce an absorption system with 
lower EW than a hotter cloud of the same column den- 
sity, which would not produce this "blanketing" . If there 
is well-developed turbulence such that at physical scales 
of the overturning regions, the column density is small in 
relation to the total column through the cloud, then this 
must be added in quadrature to the thermal 6-value in or- 
der to arrive at the total 6-value; b — y/b'f herm + b'j urb . If 
turbulence or virial motions exist on large scales, then the 
cloud should ideally be characterized by a superposition 
of independent Doppler parameters, rather than a single 
one. 

The measurement of b from the Lya line of an H I ab- 
sorber (hereafter b^ya) may reflect not only the line broad- 
ening effects of temperature and small-scale turbulence, 
but large-scale cloud motions, or physically associated sub- 
condensations; thus it is also a function of the resolution 
of the spectra. 

For clouds with no large-scale motions, the line-blanketing 
effects in clouds can be characterized by t he use of the 
relation of IChernomordik fc Ozernovl l)1993f) that approx- 
imates the equivalent width of an absorber as a function 
of optical depth r and the Doppler parameter 6; 
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The optical depth r = 1.497 x 10~ 2 \f 12 N/b, where the 
oscillator strength for Lya is /12 = 0.416, and A = 1215.67 

4 for hydrogen, b thf , rm ~ 0.128 T 1 / 2 kms' 1 
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A. Substituting, and solving for column density, 



N HI = 1.05xl0 12 6 



\ 



exp 



0.247 W 



(29) 

This relation is claimed to be accurate within 10% for all 
optical depths r < 40. For small W (W <, 150 mA), it 
reduces to the standard linear relation, 



N H i = 



W 
54.43 



x 10 13 cm -2 , 



(30) 



as noted in lPenton et~aTl l)2000b|) . 

It is possible to separate the effects of the large scale 
motions influencing the Lya line-width from the effects 
properly incorporated within the Doppler parameter by 
using observations of higher-order Lyman series absorbers 
and curve-of-growth analysis to arrive at the true b- value. 
Using the Far Ult raviolet Spectroscopic Explorer (FUSE), 
IShull let alJ <|2000|) observed 12 clouds with W > 200 mA. 
Results suggested that the intrinsic (thermal) line widths 
of the clouds were generally much smaller than that mea- 
sured from the Lya line. For the clouds with modest to 
small uncertainties in b-value, the average reduction from 
bhya to b was a factor of from 0.5 to 0.57, with disper- 
sions in the range a — 0.15 to 0.19. In studying Shull's 
results, the fractional reduction of 6-value showed no ap- 
parent correlation with cloud EW. The large dispersions 
indicate that there is a significant random element in the 
fraction of b^ ya that is due to bulk motions or clustered 
sub-condensations. 

In Fig. |S1 the derived b- values of clouds from the FUSE 
analysis are displayed as filled circles together with their 
l-<7 errors in W and b. The lower limit of their distribu- 
tion should be limited by the thermal equilibrium, as the 
turbulent broadening approaches zero. A theoretical limit 
for the thermal broadening for lower EW clouds can also 
be applied by noting the column- weighted temperatures of 
model absorption systems, as discussed in H5.2.2I and plot- 
ted as open pentagons in the figure. The apparent linear 
trend of the latter three points is extended (dashed line) 
to high EW. It is seen that the derived 6-values of large 
EW clouds obey the lower-limit that thermal broadening 
in absence of turbulence would imply. The equation of the 
line is, 

b= 11 + 0.04375 W kms -1 . (31) 

Also shown in Fig. [8] are the cloud data of lPenton et alJ 
( 2000a|), which have been reduced to 55% of their recorded 
bhya values, in anticipation of the applicability of a sim- 
ilar correction factor for lower EW clouds (see The 
appearance of a few clouds below the line of thermal broad- 
ening limit should not alarm the reader; in applying the 
mean reduction factor, we are introducing an error of or- 
der 30%; the fractional d i spersi on in the relation of b to 
bhya from the IShull etaP l)2000D data. For a cloud of true 
6-value 15 kms -1 , the 1-a error is of order 4.5 kms -1 , 
enough to displace a cloud with small bulk motions to 
well-below the line. 

Figure |H1 demonstrates that 6-values of clouds can be 
characterized as having a clear minimum for a given EW. 
Clearly, though, this will over-estimate the column density 
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Fig. 8.— The b-value data of IShull etafl l2000h using FUSE 
(large solid circles with error bars), plotted against the rest EW. 
These results imply that intrinsic 6-values for cloud s with W > 200 
mA are roughly half of the observed &Lya (see 115.31 . The cloud data 
are presented as small squares; larger tidal fields (T > 0.1) are open 
squares, and for T < 0.1, solid squares. These Doppler parameters 
have been reduced to 55% of their bL ya value, according to this 
prescription. Open pentagons represent a theoretical lower limit to 
the observed b-value, based on the column density-weighted tem- 
perature of the model clou ds for log-columns of 12.5, 13.0, 13.5, and 
14.0, left to right (see fl^Sl . The short-dashed line is an extension of 
the trend of the pentagons, showing that the FUSE data obeys the 
envelope. The long-dashed line represents the approximate median 
relationship (see text). See text for details. 



of an average cloud. However, a median b-value line may 
be drawn (long-dashed line) as a function of EW that is 
more relevant for transforming the EWDF to the CDS 
since it approximates the most probable 6-values. This 
line has the equation 



b w 17.5 + 0.042 W kms" 



(32) 



and is designed to represent the contributions of thermal 
broadening as well as bulk motions or the effects of unre- 
solved components. 

The result of applying these formulaic relations to the 
clouds via Eq. 1291 results in the mean, and void EWDFs 
(Fig. being transformed into a CDS. This is shown 
in Fig. El where the dotted lines represent the results 
of applying Eq. [21 to void, and mean cloud catalogs, 
and the solid lines represents the application of Eq. 1321 
The dashed lines accompanying these plots represent the 
CDS that would result from an untruncated isothermal 
sphere (t) = 2, resulting in a spectral slope S ~ —0.67 
top line) , and a truncated isothermal sphere (lower dashed 
line), where the truncation radius scales with r t ~ 144 h^g 1 
kpc for a cloud circular velocity v c = 21 kms -1 , about 
twice the scaled truncation radius assumed in Paper 1 
(where R t — 500X v hj^ kpc). At this distance the over- 
density of baryons and DM (assuming baryons follow DM) 
is 2.57 fib Pent, and thus shows the strong effect of a trun- 
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cation (such as from ram pressure stripping of diffuse, dis- 
sipated gas in filaments) on the CDS slope at low N HI . 
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F ig. 9. — The observed c olum n density spectrum derived using Eq. 
1311 (dotted lines), and Eq. 1321 (solid lines), for the volume- weighted 
column density distribution (upper pair of lines) and the void CDS 
(T < 0.03, lower pair). These distributions are compared to an 
analytic fit of arbitrary normalization for an untruncated isothermal 
baryon distribution (upper dashed line) and a truncated isothermal 
model with rt ~ 144 kpc for a halo with v c ~ 21 kms — 1 (lower 
dashed line). The effect of truncation is to diminish the number of 
low column density absorbers. The slope of the untruncated line is 
S = -0.67. 

For the purposes of transforming an EWDF into a CDS, 
or vice versa the median relationship of 6(W) (Eq. I32[l will 
henceforth be used. Though the error at column densities 
Nhi J> 10 13 ' 5 cm~ 2 may be substantial, smaller columns 
will be minimally affected, as can be seen in Fig. [5] 

6. MODELING DISTRIBUTIONS OF SUBGALACTIC HALOS 

The foregoing developments now enable the modeling of 
the CDS with distributions of model clouds. Cloud profiles 
can be combined into a column density spectrum, given a 
distribution of halo velocities. The primary method of de- 
termining the halo distribution function used here involves 
transforming the luminosity function (LF) into a halo ve- 
locity distribution function (hereafter HVDF). Additional 
useful information can be garnered by comparing our re- 
sults to those emerging from numerical simulations. 

The required procedure can be outlined as follows: First, 
extract the column density of a model cloud as a function 
of projected radius from az = product of the 1-D simu- 
lation. Secondly, combine this with a distribution of halos, 
and compare this synthetic spectrum to observations. 

Since clouds are here assumed to be sub-galactic pertur- 
bations, I approach the subject by transforming the LF of 
galaxies into a HVDF that can be used to synthesize a 
CDS out of model clouds. The halo mass and velocity dis- 
tribution functions of galaxies are determined by applying 
the Tully-Fisher (TF) relation to the galaxy luminosity 
function. 



The TF relation be used to transform the luminosity 
function into an HVDF. It can be represented in many 
ways. Perhaps the simplest is as the power-law, 



C_ 







(33) 



where v* = 161 krns" 1 (see §5.1 of Paper 1). The value of 
the constant is a function of the photometric band being 
used. For the .B-band, = 2. 91, and for the if-band, 
= 3.51 l)Tullv fc PierceL 1200(1) . Note that all symbols 
with a star, e.g., M* , R^, & etc., refer to expectation 
values for an C* system, itself referring to the point of 
inflection of the Schechter form of the LF (i.e., the "knee" ) . 

As long as there is a calibrated TF relation correspond- 
ing to (i.e., in the same spectral bandpass as) a given LF, 
one may apply a transformation of the LF into a mass, 
or velocity-DF. Below, is detailed the application of this 
transformation. 

6.1. Transforming the LF to a Velocity or Mass DF 

The halo velocity or mass distributions are based on 
the luminosity function. For these distribution functions I 
make the following syntactical simplifications, 



X, 



v c 

x v = —, 
v* 



(34) 
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M 

M* 



where M. represents the halo mass. The Schechter lumi- 
nosity function is, 



(35) 

l)Schechteiill976j) and is a differential DF. A 1-D analogue 
to the Jacobian will transform the LF into something else. 
For instance, to transform a Schechter distribution </>£ into 
a HVDF (cj> v ) we have, (f> v — 4>cdXc/dX Vc . Partial deriva- 
tives are calculated from the following scaling relations for 
a truncated ISOT halo (see Paper 1, Eq. 29), 



X„ 



Rt 
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(36) 



where Rt is the galaxy halo truncation radius. The trun- 
cation radius utilized here, for clouds as well as galaxies, is 
such that the total mass density at R t is p(Rt) = 10 p cr it- 
For the NFW halo the following scaling relations obtain 
(Paper 1, Eq. 40), 



X v — 



Rn 
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'-'ft 



0.67 
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V/3 



_ Y 0.30 



(37) 



As with Paper 1, we correct both the ISOT and NFW 
halos for the amount of baryons expected to be associated 
with it according to the assumed value of Qb/Sl m . With 
baryons included, v* — 161.5 kms -1 (see §5.1.1 of Paper 
1), but the dark halo has a halo v* = 144 kms -1 . 
Since X c = X$ (Eq. E3, 



dXc 
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This halo velocity distribution function is identical for 
ISOT and NFW halos. The mass function for the ISOT 
halo is, 

, , dX c (P4>*c\ „/3( a +l)/3-l -x 13 ' 3 tm\ 

while for the NFW halo we have, 

, , 9X C , „,» . v 0.SP(a+l)-l -X°f e fAn \ 

Note that the exponential may be safely ignored for clouds, 
that have X v and Xm -C 1. 

The normalizations for each of these distributions can 
easily be calculated using the normalizations of the lumi- 
nosity function. For the B-band LF, it is thought wise 
to interpolate the bracketing Sloan Digital Sky Survey 
(SDSS) normalizations to the B-band owing to the ap- 
parent consistency of their normalizations (see Appendix 
of Paper 1) and the large database used. We therefore use, 

cj)* B = 0.022 h 3 Mpc -3 = 0.0093 h^Mpc" 3 . (41) 

The distribution functions given in this section provide 
the basic framework for placing our halos into the universe. 
Although there are some conflicting indicators of which 
value of the slope a to apply, we may either prevail upon 
the criterion that a should be steep enough to explain the 
void cloud line density, or find a preferred value for a based 
on other compelling work. 

6.2. Comparison with Previous HVDF Simulations 

Cosmological sim u lations of hiera rchical clustering (e.g., 
iKlvnin et all Il999t iKlvpinl |2002|) (hereafter K99, and 
K02, respectively) provide the basis for deriving veloc- 
ity and mass distribution functions of halos. How do 
the slopes of these distributions relate to the functional 
forms derived above? K02 noted an environmental depen- 
dency in the halo velocity distribution, with slopes —3.1 in 
groups, and —4.0 for isolated halos. With these facts we 
can find the Schechter faint-end slope parameter a that 
is consistent with these values. For small halo velocities, 
the HVDF (Eq. 1381) can be simplified by removing the ex- 
ponential which is ~ 1 for cloud halos. The HVDF is the 
same for ISOT and NFW, but there are two different envi- 
ronments; the dense (i.e., group), and the isolated environ- 
ments. I alternately equate the HVDF slope, f3(a + 1) — 1, 
to —3.1, and —4.0, and derive the slope parameter us- 
ing Pb — 2.91: a = —1.72, and —2.03, for group, and 
isolated halo distributions, respectively. Note that these 
simulations presume NFW halos. It stands to reason that 
the isolated halo distribution would be more applicable to 
void cloud distributions. 

By equating the mass function slopes in Eqs. 1391 and I4TH 
(excised of their exponentials) to that of the differential 
mass function with a slope of —1.8, as found by K02, 1 find, 
a = -1.83 and -1.93, for ISOT, and NFW, respectively 
(/3b = 2.91). These represent values a mixture of isolated 
and group halos. 

If isolated NFW halos have a slope parameter a = 
—2.03, perhaps the slope slope parameter for isolated 
ISOT halos can be estimated by the difference between 



ISOT and NFW slopes for mixed environments (a = —1.83 
and —1.93, respectively, as above). Since the ISOT slope 
was flatter than the NFW mass function by 0.1 for mixed 
environments, then this would imply that for isolated 
ISOT halos, a = -2.03 + 0.1 ~ -1.93. I consider a slightly 
rounded value a — —1.95 to be a provisional order-of- 
magnitude prediction for ISOT clouds in voids. 

This steep slope, in co mparison to a faint-e nd slope 
of a w —1.25 for galaxies ijBlanton et all l2001jh demon- 
strates the source of the "missing galactic satellites" prob- 
lem (e.g., K99). The suggestion of K99 is that the missing 
satellites may be clouds, rather than dwarf galaxies. In 
this case, their work may imply that the parameter a, if 
it is taken as referring to the demographics of halos in 
general, rather than the luminosity of halos, may be in- 
dicative of the presence of numerous small halos that are 
not luminous. Whether these are LSB galaxies, HVCs or 
Lya clouds, a steep EWDF slope parameter suggests that 
a significant amount of mass may be tied up either in small 
halos, or in an unconsolidated component, as we shall see 

in ina 

6.3. Assembling the CDS 

The column density spectrum is calculated from model 
clouds according to the equation, 

= f f «r% H1 6X v (n), (42) 

a z -Ho — ?, 

n— 

where 4>(v c ) is given by Eq. |2H1 r N HI is the radius at which 
a column density Njji is detected, X v (n) is given by Eq. 
1281 and 8X v (n) is the bin width. The cadence in X v indi- 
cates that SX V = V„(10 025 -10~ 025 ) = 0.1152 X v . The 
radius at which a given column density can be observed 
in a cloud of a given v c is determined by the method de- 
scribed in H5.2.1I The model column density spectra can 
then be compared with the observations. 

The goal of modeling the distribution of clouds is to re- 
produce the observations. Not only must the line densities 
as a function of column density be matched, but it is also 
important to reproduce the shape of the CDS. 

Equation |2H1 shows how the halo distribution function is 
calculated as a function of X v , and demonstrates the role 
of the faint end slope a and the TF slope (3 in determining 
the low-velocity end of the HVDF. Since it is theorized 
that the distribution of sub-galactic halos is in continu- 
ity with the galaxy halo EWDF, the normalization of the 
galaxy luminosity function in voids will be a good place 
to start the modeling; the modeled distribution of halos 
must be based on a normalization <p*(void). It has been 
suggested that galaxies are present in voids at ~ 1/10 the 
average density l)El-Ad fc PiranL I1997LI2000D. A somewhat 
larger galaxy densit y is found in Bootes l|Dev et all Il995t 
ICruzen et all 12002^ . An upper-limit of current thin king 
might be represented by the simulations of lCen fe Ostrikerl 
(1999), which assumes that fl m — 0.4 and estimates that 
f2y/f2 m ~ 0.26, consistent, perhaps, with roughly an 
equal mass in clouds and in galaxies. A low normaliza- 
tion would seem to imply a low total void density, yet 
since the spectral slope of the void cloud HVDF may be 
steep, this is not necessarily inconsistent with a large den- 
sity. Therefore, I provisionally adopt the normalization, 



C. V. Manning 



15 



cj)*(void) = 0.10* = 9.3 x lO^h^Mpc 3 (see Eq. HI) , 
with hopes that it can explain the data. 

The somewhat arduous transformation of the TF re- 
lation to an HVDF, and the number of parameters that 
need to be evaluated, may give the impression that there 
are many free parameters utilized in this analysis. This 
is a misconception, for none of the parameters are given 
arbitrary values. The slope parameter a is reasonably 
constrained by the modified Press-Schechter simulations 
of K02 (see ^6.2[1 . the TF slope (3 is constrained by 
the B-band TF relati on, which has a well-known value 
(Tul lv fc Piercel IToOff). and the normalization of the void 
EWDF is observationally constrained, as noted above. Of 
these, the latter would seem to have the greatest range of 
uncertainty, and is the only parameter used as though it 
were a truly free parameter. However, none of the values 
of these parameters are arbitrary, or disconnected from 
observations. 

The effects of using various slope parameters a on the 
derived CDS are discussed in the following section. 

6.4. The model CDS with slope parameter a 

A constant slope parameter a is the most obvious be- 
ginning point for investigating the relationships between 
our parameters and the product CDS. 

6.4.1. The ISOT model 

I first consider the ISOT model. Recall that with the 
ISOT halo, clouds with large v c significantly restrain the 
expansion of baryons from the dark halo f US.l.lf) . As a 
result, columns Nhi > 10 124 cm -2 are detected for the 
velocity range 10 < v c < 25 kms -1 . Fi gure 1101 displays 
model CDS (solid curves) for various slope parameters a 
using the above normalization and the ISOT model profiles 
(Fig. [SJ). When applying slopes of a = —1.0 or —1.4, 
low line densities result. If the slope is —1.4, then the 
normalization will have to be increased by a factor of 30 
in order to give rough agreement with the observed line 
density, making the normalization in void space greater 
than the average by a factor of 3, which is absurd. When 
a w —1.95, as suggested in H6.2I the resulting modeled 
line density is comparable to the observed value. In order 
to fit the low column density end of the CDS, a slope 
parameter of order a ~ —3.5 is required (dotted line in 
Fig. 1101) . However, integrating over the HVDF (for details, 
see H6.6f) . the total mass density for clouds 5 < v c < 25 
kms -1 is il c i ;> 0.028, of order 10% of the assumed mean 
density Q m = 0.3, when the truncation radius is taken to 
be r t = 500Xu h^ 1 kpc, where X v is given by Eq. 123 as 
was the standard with Paper 1 (§5.1.1). When systems of 
velocity 10 < v c < 150 kms" 1 are used, then fly = 0.09, 
and fV/^m — 0.3. The truncation radius is such that the 
density at r t is about 10 times the critical. Considering 
that the rest of void space is unlikely to be absolutely 
empty, this value for Qy must be a lower limit. 

6.4.2. The NFW model 

NFW halos do not produce any absorbers with Nhi > 
10 12 - 4 cm" 2 in the range 10 < v c < 22 A kms -1 ; the next 
higher velocity bin (n = 13 in Eq. 20) experiences col- 
lapse in the center, presumably forming stars. The line 
density, even at Nhi = 10 122 cm -2 is less than one cloud 
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Fig. 10. — A comparison of the observed void column density 
spectrum for low-2 clouds (T < 0.1) (short-dashed line) with spectra 
that would result from a HVDF resulting from an extension of the 
SLF faint-end slope a = —1.0, —1.4, and —1.95, normalized at 
0.10* (solid lines). To fit the slope at Nhi < 10 13 3 cm" 2 , an 
HVDF with a slope parameter a ~ —3.5 is required (dotted line), 
though the required normalization is an order of magnitude lower 
than 0.10*. 

per unit redshift for parameters that produce an accept- 
able fit with ISOT halos. The reasons for this shortfall are 
two-fold. First, NFW halos are about 5 times less massive 
than an ISOT halo of the same v c , and hence, have 5 times 
fewer baryons. Second, the lower mass of the halos can- 
not as effectively slow the baryonic outflow and dispersal 
as can the ISOT model; the trend of the restraint of the 
outflowing gas for large v c , noted in the inset of Fig. H3 
is not apparent in the NFW profiles (Fig. [7). Further, 
comparison of Figs HJI and [7\ shows that the central bary- 
onic densities of NFW halos are of order 1/10 of an ISOT 
halo at comparable halo velocities. Since no absorption 
systems are produced, this halo model fails to explain the 
observed EWDF, and I forego any further consideration of 
it. 

6.5. The Shape of the CDS 

In this section I take a critical look at the observed CDS 
to see what degree of confidence is justified in its slope at 
the low, or "high" column density end (~ 10 14 cm -2 ). 
Figure ^| shows the Poisson error bars for data at z < 
0.036 (see §6.2.1 of Paper 1). Adjacent error bars are 
not independent of each other since the same clouds may 
contribute to a range of EWs. When data from all red- 
shifts are used, the Poisson errors are smaller, but the 
slope at the low column-density end (W <J 40 mA, or 
Nhi <j lO 13 3 cm -2 ) is moderately flatter than if only 
data from z < 0.036 is used (see Fig. 11 and 12 of Pa- 
per 1). This is probably attributable to an underestima- 
tion of the tidal field at higher redshifts, with consequent 
contamination of void catalogs with non-void clouds, all 



16 



Modeling the Void H I Column Density Spectrum 



occasioned by the probable non-detection of less-luminous 
galaxies in the proximity of the clouds. I conclude that 
the steepness of the slope apparent at log(iVffj) <; 13.3 in 
Fig. I1UI is highly significant, though what it means is yet 
to be determined (see ^6.611 . 

However, the "bump" in the CDS at columns Nhi <; 
lO 13 5 cm -2 , remains essentially unchanged when higher 
redshift data are excluded. Is this a real feature? There 
are five clouds with Njji > 100 mA (which, using Eqs. B2l 
andl2*9l implies N H i <; 10 1335 cm -2 ) contributing to this 
end of the low-z void distribution (T < 0.1), as opposed to 
n c i ~ 20 at Nhi < 10 12 9 cm -2 . It is helpful, at this point, 
to take a detailed look at the five clouds with W > 100 
mA in the low-z data with the help of a program that as- 
sesses the prospects for the physical association of clouds 
(i.e., gravitationally bound) with galaxies on the basis of 
the galaxy absolute magnitude (mass is derived assuming 
the more massive ISOT halo), the relative velocity of cloud 
and galaxy, and the projected radius, (Manning, in prepa- 
ration). The galaxies considered for this routine are those 
within 750 kins -1 of the cloud in question if the projected 
radius r p < 2 Mpc, or within 350 kms -1 if outside this. 

Two of the five low-z clouds have derived local tidal 
fields of T = 0.073 and 0.098. Both are found to have one 
sub-£* galaxy at a projected distance r p > 2 Mpc, with 
cloud-galaxy LOS relative velocities of 66 and 72 kms' 1 , 
respectively, for which the program predicts probabilities 
of physical association of ~ 13%, and 4%, respectively. 
These clouds can be found in Fig. 4 of Paper 1 in the sight- 
lines to MRK 335 (Fig. 4k, 1241.09 A) and AKN 120 (Fig. 
4h, 1247.94 A), respectively. But even this possible phys- 
ical association could not produce a significant peculiar 
velocity along the LOS because of the large projected dis- 
tance to the galaxy, and hence the attributed tide cannot 
be far wrong. Of the other three, two have tidal fields T w 
1 x 10 -3 , and have no galaxies within the search area. The 
other, with T ~ 7 x 10 -3 , has a late-type sub-£* galaxy at 
a pr ojected radius r„ ~ 430 kpc (the CfA Redshift Cat- 
alog, fauchra et al.lll990L iHuchra. Geller. fc Corwinlll995l 
i Marzke. Huchra. fe Gellerll996LlGrogin. Geller. fe Huchral 
Il998l IHuchra. Voeelev. fc Gellerlll999|) 5 . but at a relative 
velocity of Aw ~ 350 kms -1 . The probability of associ- 
ation was rated at ~ 4%. The galaxy velocity error was 
listed as Sv = 1 kms -1 , and the cloud redshift error from 
iPenton et al.l l)2000af) was listed as Sv = 6 kms -1 . Thus, 
these five clouds have few prospects for a physical asso- 
ciation in which significant peculiar velocities from grav- 
itational interactions could result, and fewer chances for 
associations that would give significant LOS peculiar ve- 
locities. With high confidence, I claim these 5 clouds are 
truly remote, and are not anomalously present in void cat- 
alogs. 

It is noteworthy that the void filling factor at T ^ 0.01 is 
« 0.44 (see Fig. 7, Paper 1), while at T < 0.1, f v ~ 0.86. 
Thus the volume in which the three low T clouds appear 
is almost exactly the same as the volume in which the two 
at 0.01 < T < 0.1 clouds appear, making it seem likely 
that large clouds may be more prevalent in deep voids 
than in the void periphery. But the statistics are small, 
and we would be wise to seek more data before forming a 

sec http://cfa- www. harvard. edu~huchra/ 



conclusion. 

A possible explanation for those clouds that produce 
the excess above the model at Nhi ~ 10 13 7 (see Fig. 
Ill) [I is that there exist clouds/proto-galaxies in voids with 
halo velocities v c > 25 kms -1 , for the analysis of the 
CDS of individual clouds shows that the larger the circular 
velocity of a cloud, the flatter the slope of the CDS that 
it produces; such additions would produce more absorbers 
at log Nhi ~ 13 to 14 ( cm -2 ). This can be understood by 
noting that the more massive the cloud, the more strongly 
it will hold its baryons, producing a steep density gradient 
at the edge, and resulting in a flat CDS (c/. Eqs. [TBI and 
[T7|. 

Recall that the gas in larger halos (v c > 25 kms" 1 ) col- 
lapse and produce stars in the 1-D simulations. However, if 
clouds have grown from hierarchical merging during their 
evolution, they would have been smaller during the early, 
high-density stages of evolution, achieving their larger size 
only recently. Because of the early evaporation of baryons 
from their core region, there might not be sufficient bary- 
onic mass in the core to produce the cooling needed to 
trigger star formation. Therefore, with growing halos, the 
limiting circular velocity (as measured at z = 0) for SF is 
probably higher than 25 kms -1 . 

But before pursuing this line of thought, it is impor- 
tant to consider the repercussions of taking literally the 
significant steep CDS slope at Nhi <> 10 13 ' 3 cm -2 noted 
above. 

6.6. Total mass distribution and slope parameter a 

As noted in ^6.4.11 the low column density end of the 
void CDS can be well-fit by an HVDF with a very steep 
slope parameter a w —3.5. In this section I inspect the 
physical repercussions of this. As argued in H6.21 the 
HVDF may be simplified by removing the exponential (Eq. 
1351) when X v <C 1. Consider logarithmic bins of v c such 
as was used in the simulations here; construct the bins 
such that each bin has a width 5 X v = X v Svq/vq, where 
3vq/vo <C 1 is some appropriate constant. Then from Eq. 
1381 the number density of systems within a bin is, 

dn U n = KX^ +1) - 1} xJ^\. (43) 

Thus, the comoving number density per bin scales as 

Sn bln ocX^ a+1 \ (44) 

When a = — 1, each logarithmic bin has an equal number 
of systems in it. However, when we consider the mass 
in each bin, we must multiply the number density by the 
mass per object. 

For the ISOT halo, the mass scales as M oc V 3 (Eq. 

m, 

SM bm oc AT^ +1 )+ 3 oc X^ a+1)/3+ \ (45) 

that has equal masses in logarithmic bins when a = — (1 + 
3//3) ~ —2.03 for /3g = 2.91. For this value of a, the mass 
function is flat in logarithmic bins. 

For the NFW halo, virial mass M 20 o oc v 3 „ 3 ° Eq. I37|) . 
The NFW bin mass is a somewhat different relation, 

SM bm oc Xf)(« +1 )+ 3 - 3 , (46) 
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that has equal masses in logarithmic bins when a = 
— 3.3//? — 1 ~ —2.13 for the same fl B . 

Recall that a CDS formed using ISOT cloud profiles and 
a HVDF with a = —3.5 would have an increasing mass in 
bins of progressively smaller velocity v c , implying a diver- 
gence of the integrated mass. The sobering implications of 
this fact are softened, however, if sub-clustering is occur- 
ring. Further, it is reasonable that sub-halos should form 
in the flanks of larger halos, for the time-scale for collapse 
of perturbations is proportional to the inverse square root 
of the density, making it more likely for them to form 
there than in the regions between clouds. Thus, many of 
the small clouds that are detected in the CDS may not be 
independent, isolated clouds, but loosely associated with 
larger clouds. Since they are fashioned out of the stuff of 
the main cloud, they are not taken to add to its mass, un- 
less such sub-halos are found to lie outside the truncation 
radius Rt- If sub-condensations exist, then accretion or 
hierarchical agglomeration is a distinct possibility in the 
halos of H I clouds in voids. 

7. EXTENDING MODEL CLOUDS TO HIGHER V c 

Up to this point, the halos that have been used were 
fixed at z = 6.5. Such halos remained star-free only up 
to circular velocities of only up to v c w 25 kms -1 (SJSJ. 
However, the steep spectral slope of the void CDS at low 
column density appeared to require the presence of sub- 
clustered halos, since otherwise the indicative mass func- 
tion would diverge. Lik ewise, it is difficult to reconcile 
the existence of CHVCs llBraun fc Burtonl 120001) w ith H 
I mass approaching 10 7 (|de Heii et all l2002|) . with 
model clouds having maximum H I masses only <; 10 6 M Q . 

The behavior of the observed CDS relative to model 
CDS thus suggests the possibility of continued accretion 
and halo growth during the period from z = 6.5 to the 
present; that sub-halos cluster about larger halos implies 
that the circular velocities of clouds could grow by the 
accretion of those halos; their growth would explain the 
presence of CHVCs 6 with H I masses approaching 10 7 M@. 
Figure ITT1 shows H I mass for fixed halos (open circles), 
corresponding to each circular velocity bin) . The trend of 
fixed halos shown in Fig. 111! halos of velocity v c ;> 40 
kins -1 could reach this H I mass. 

Recall, however, that attempts to run models of such 
large halo velocities failed for fixed halos, because the cen- 
tral baryonic core, thought to roughly follow the dark mat- 
ter before reionization, was not able to support itself, as 
described in However, if halos are slowly "grown" , such 
as by accretion, then larger halo velocities at z — may 
be attainable without these indications of star formation, 
with H I masses more in line with observations of CHVCs 
may be attained. 

In what follows, I adopt this model of accretion. A rough 
approximation of that growth rate was made by assuming 
that the largest fixed cloud that escapes collapse into star 
formation, could grow to the largest cl oud which today ha s 
an H I mass of 10 7 M Q , as asserted bv lBraun et al.1 l)2002j) . 
The circular velocity of that halo was approximated by 
connecting the trend in H I mass with fixed-halo circular 

6 the optical surface brightness of local CHVCs is constrained to 
be < 10% of the faintest Local Group dSph (L. Blitz, private com- 
munication 2002) 



60 

3 5 



o 'Fixed' halos 


i i 

- 


• 'Grown' halos 
Js 


- 
- 


0/ 


- 


of 




J 




- 





0.8 1 1.2 1.4 1.6 

Log(v c ) (km/s) 

Fig. 11. — The H I mass in te rms of the circular velocity for 
clouds that grow according to Eq. 1471 (filled circles and solid line) 
compared to that of the fixed halo (open circles, dotted line). The 
logarithm of the H I mass is shown as a function of the logarithm 
of the cloud circular velocity (solid line) for untruncated clouds. 



velocity to the maximum cloud H I mass of 10 7 M (see 
Fig. Ell- 




Fig. 12. — Density profiles for clouds with circular velocity v c ~ 
9.4, 15.0, and 23.7 kms" 1 for fixe d halos (dashed lines) and for halos 
which grow according to Eq. 1471 (solid lines). The mor e power- law 
like profile results in a steeper CDS slope (see Fig. 1131 . 

The clouds are "grown" gradually, that is, by inserting 
increments of mass into the Lagrangian bins at each time 
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step in a way that mimics the equation, 

v c {t)=v c (0)z ka \ (47) 

where v c (Q) is the circular velocity at z = 6.5, the accre- 
tion constant k a « 0.78 in units of inverse Hubble times 
tjj , and t is in units of tn- As DM is gradually loaded 
into the DM bins, their physical extents are scaled-up to 
mimic the larger core-radius of larger clouds (see §5.1.1 of 
Paper 1 for scaling relations). Mass is added to the bary- 
onic bins in the same way, however, the radial positions 
are not altered directly, but allowed to evolve according to 
the pressure gradients and total gravity of the cloud. This 
approximation may result in a slightly more concentrated 
distribution of baryons than if the gas was originally as- 
sociated with the accreting sub-halos. But the inward mi- 
gration of DM would tend to pull the gas inward anyway. 
It appears that the only significant difference between the 
model and the envisioned reality may be a weaker discon- 
tinuity in the model between the outflowing wave from the 
original ionization of the cloud and the infalling gas from 
the periphery. But the relatively mild difference between 
fixed halo, and grown halo evolution suggests that a thor- 
oughly accurate treatment of this problem would present 
a final profile much like the grown halos. 

Results for some of the halos are shown in Fig. ^1 where 
the solid density profiles represent the "grown" halos, and 
the accompanying dashed lines represent fixed halos of the 
same circular velocity. In general, the resulting halos have 
a modestly lower central density and a modestly greater 
density outside the core. In this way, clouds with central 
bin velocities up to v c = 33.5 kms -1 were generated with 
no central baryon collapse, three more velocity bins than 
with fixed halos. For still larger halos, however, the inner- 
most bins collapse, indicating the plausible formation of 
stars. The bin with center 37.6 kms -1 collapsed at very 
low redshift, indicating the possibility of very young void 
galaxies forming from such clouds. Thus, the growth of 
void clouds could expl ain the existence of blu e compact 
galaxies at voids edges l(Pustil'nik et all Il995lk with halo 
circular velocities v c J> 35 kms -1 . 

For a given normalization <fi v , extending the halos to 
larger circular velocities in this way results in a larger line 
density at a given column density. As a result, the best fit 
to the observed CDS requires lowering the normalization 
to <py ~ 0.06 4>* . The resulting CDS, using the grown ha- 
los with HVDF slope parameter a = —1.95, is presented 
in Fig. ^|as the solid line; the dashed line represents the 
fixed halos with cf> v = 0.1 cf>* . Note that the quality of 
the fit is greatly improved at large column densities. The 
excess at small column density is still expected, due to the 
sub-clustering of small halos in the flanks of major clouds. 
Modeling this component requires semi-analytic work be- 
yond the scope of the present paper, and is deferred. 

If the mass function is integrated over the range of 
model halo velocities which are detectable, the mass den- 
sity in voids is, Qy > 0.022, for <\>* v = 0.06 0* and 
R t = 500 X v kpc. This is 78% of the mass density for fixed 
halos. The reason for this difference is first, the lower nor- 
malization, but this is partially compensated by the larger 
halos which emerge from the grown halo models. 

8. DOPPLER PARAMETERS OF MODEL ABSORBERS 




12.5 13 13.5 14 

Log(N HI ) 



Fig. 13. — Predicted CDS using an HVDF with slope param- 
eter a = —1.95 for clouds evolved with fixed halo s (d ashed line), 
compared with halos which grow according to Eq. 1471 The latter 
halos are detectable in 4 more bins than the former, one centered 
at 9.4 km/s, and three at 26.6, 29.9 and 33.5 km/s. The normal- 
ization for the grown clouds was (fry = 0.06 <fr* , and for fixed halos, 
^=0.10*. 

Recall that the diffuse absorbers expected from N-body 
simulations predict large Doppler parameters, and that 
this was one of the reasons that this author abandoned the 
diffuse cloud model (JQ|. The evolved baryon distributions 
pertinent to a given halo model have a wealth of detail 
in the program output of the 1-D simulations. This sec- 
tion investigates the Doppler parameters of the model ab- 
sorbers. Fig. 1141 shows the cumulative line density dAf/dz 
contributed by each of the model halos, as normalized by 
0.1^* for the fixed halos (panel a), and 0.060* for grown 
halos (panel b); both with slope parameter a = —1.95. 
Note that low column density absorbers are produced in a 
wide range of cloud models. 

Figure ^| shows histograms of void (dotted) and non- 
void (solid lines) clouds for W < 50 mA (panel o) and W > 
50 mA (panel 6), for the low-redshift sample, z < 0.036. 
Panel a shows that low EW void clouds have remarkably 
uniform 6-values, with a near-symmetrical profile with a 
maximum at ~ 30 kms -1 , while panel b shows that of the 
higher EW clouds, non-void clouds predominate only at 
higher b-values. That is, almost half of the clouds with 
W > 50 mA and v c > 40 kms -1 are void clouds. 

Fig. I14t > shows that, for grown clouds, a cloud of column 
density Nhi = 10 124 cm -2 can be produced by clouds of 
circular velocity 10.6 < v c <; 35 kms -1 . For purposes of 
illustration, I now undertake to closely inspect the thermal 
and velocity structure of model absorption lines produced 
by the extremes of halos producing this column density. 

I estimate the core temperature of clouds with halo ve- 
locities v c = 10.6 and 33.5 kms -1 using the method of 
^5.2.21 The smaller cloud is found to have a weighted 
temperature of 7863 K, while the larger cloud has a tern- 
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Fig. 14. — Panels a and b show the cumulative line density 
(greater than a given column density) as a function of the halo 
velocity bin for fixed, and grown halos, respectively. The top line is 
the minimum detected column density, Njjj = 10 12 ' 4 cm -2 , with 
subsequent lines representing columns larger by 0.2 in log Nhi- A 
few of the lines are labeled by the log of their column density. Low 
column density clouds are produced by clouds with a wide range of 
halo velocities. 
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Fig. 15. — Cloud Doppler parameter histograms from 
IPenton et ahl 12000afl GHRS data for void (dotted) and non-void 
(solid line) catalogs (using z < 0.036), for clouds with W < 50 mA 
(left-hand plot), and W > 50 mA (right). For low-EW clouds there 
is a narrow range of b- values, while for larger clouds, b- values ap- 
pears to vary over a wide range. Notice the difference in vertical 
scales between panels a and b. 



perature of 5375 K. This can be understood by reference 
to Fig. for the larger cloud has a much larger impact pa- 
rameter than the smaller cloud (r p ~ 500 kpc versus ~ 55 
kpc, respectively), and hence has lower average densities. 

There is apt to be large-scale turbulence in the clouds 
following reionization, especially toward the center where 
the density, and its gradient remain high, and the potential 
for convective instability is large. A plausible upper-limit 
for a cloud's turbulent velocity is the cloud virial velocity. 
The virial temperature is, 



T 



1 vlnrriH 

2 k b 



(48) 



and may be converted to a virial &-value by b V i r = 
V0.0165T„i,.. The virial temperatures of the two clouds 
(T vir = 6776 K, and 67,756 K, for the smaller, and larger 
clouds, respectively) may contribute to the Doppler pa- 
rameter. For the small 10.6 kms -1 cloud, the addition to 
b contributed by turbulence is apt to be small since the T V i r 
is less than the thermal temperature. For the large cloud, 
the turbulent component could dominate, depending on 
the impact parameter (hence column density). Since the 



impact parameter of the large Nhi = 10 1 cm" 2 cloud 
is large, the sightline won't penetrate the virializing por- 
tion of the cloud (r p <^ 125 kpc). Instead, the velocity 
field of the expanding cloud may be an important addi- 
tion. The LOS velocity difference within the FWHM of 
the line is found to be 25 kms" 1 over the LOS distance, 
h/2 — 534 kpc, yielding b ve i ~ 15 kms" 1 . If this material 
were in the Hubble flow, the FWHM velocity would have 
been ~ 40 kms -1 , showing the physical effect of the neg- 
ative peculiar velocities previously noted for large clouds 
(see Fig. gj. 

The 6-value associated with the smalle r Nhi = 10 12 4 
cm" 2 cloud is then b = \/b'f he 



~bj urb < 15.5 kms" 



where the b ve i component is ignored since it is found to 
be less than 2 kms -1 within the FWHM of the line. The 
direction of the inequality in the equation stems from the 
fact that v t urb is an upper limit for bulk motions in these 
absorbers. The lower limit is the thermal 6-value; b = 
11.4 kms" 1 . For the larger cloud, b = \Jb'* herm + b 2 vel > 
17.7kms _1 , where the turbulent (virial) motions of the 
cloud have not been included since only the outer cloud 
is probed. This velocity is a lower limit because some 
turbulent motions are likely. 

Larger clouds, such as Nhi = 10 13 5 cm" 2 (W ~ 130 
mA), have averaged temperatures on order ~ 16,900 K 
(btherm = 16.7 kms" 1 ), produced by halos of velocity 
17 < v c < 33 kms" 1 . They give total fe-values b <J 24 
to b <, 37 kms -1 , where I have ignored the cloud ex- 
pansion broadening, which is dwarfed by possible turbu- 
lent motions. The impact parameters for these systems, 
120 <; r p <J 180 hyg 1 kpc, may probe virializing regions. 

These v alues agree well with the range of b- values pre- 
sented in l|Dave &: Trippl l2001[) (see their Fig. 6a). 

In c omparing the analysis of b- values in iDave fc Trippl 
l)2001|) to those in this paper, it is interesting to note 
that though the former predict higher temperatures for 
a given gas density (see Fig. |SJ, the thermal b- values as 
a function of a given column density (Eq. |SJ) are smaller. 
Table 1 shows a compilation of thermal 6-values; column 
1 is the log column density; the next column shows the 
IDave &: Trippl l)2001j) thermal 6- values; third are those of 
the present study. It would appear that at the same col- 
umn density, the former must have substantially lower H I 
densities (by a factor of ~ 3 to 5) than the discrete clouds 
considered here, and therefore have proportionately larger 
integration paths. The irony sh ould not be missed; the 
absorbers of lDave fc Trippl l)2001|) are supposed to occupy 
a dense environment, but are less dense than the void ab- 
sorbers modeled in the present effort. 
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I have shown that centrally condensed clouds based 
on sub-galactic halos produce Doppler p arameter distri- 
butio ns consistent with observed clouds l|Dave fc Trippl 
1200 lft . Per haps, however good agreement is not ex- 
pected since iDave fc Trippl l)2001j) purport to study pre- 
dominantly non-void clouds, while the present study mod- 
els void clouds. On the other hand, if the lines of s i ght coy - 
ered in their paper are like those in lPenton et alJ l)2000alK 
then Fig. 1141 would indicate that most clouds, especially 
those of lower column density, are in fact void clouds, and 
substantial agreement with the echelle observations would 
be expected. 

However, it is unlikely that one could decide between 
giant sheets and discrete, sub-galactic absorbers on the 
basis of Doppler parameters alone. The conclusion to be 
taken away from this section is that discrete halos produce 
results that are consistent with observations of resolved 
Doppler parameters in low-redshift Lya clouds. 

9. SUMMARY 

In the introduction, two possibly problematic aspects 
of recent N-body/hydro simulations were noted; first, the 
fact that Lya absorbers are expected to be produced in 
diffuse clouds by a "fluctuating Gunn Peterson" effect 
(Croft et all 119981 IDave fc Trippl |2001(K rather than by 
absorption in discrete halos. Second, that low-redshift ab- 
sorbers are expected in diffuse filamentary structures sur- 
rounding galaxy concentrations, but not in voids. 

The analysis in |2| showed that indeed, diffuse structures 
could not provide an explanation for the void absorbers. 
The analysis in the balance of the paper showed that dis- 
crete halos may restrain the baryons from evaporating to 
a degree sufficient for absorption systems to be produced 
that are consistent with observations. This implies that 
it is not a good approximation to assume that small ha- 
los would not retain a significant quantity of their baryons 
following reionization. 

However, the problem goes farther than just contest- 
ing the N-body Lya cloud simulations. For the current 
halo paradigm, NFW, also fails to explain the absorbers 
since those halos small enough to escape star formation 
are not massive enough to hold baryons against photo- 
evaporation. In short, a more massive halo is needed. 

I have shown that an isothermal halo can explain the ab- 
sorbers, within the cosmological parameters investigated 
here. Some might argue that the ISOT halo is not physi- 
cally realistic since they cannot be virialized. However, I 
am not proposing that "ISOT halos" are totally virialized. 
Instead, I propose (as noted in H4. 1.4(1 that the isothermal 
density profile is a fair description only of the total mass 
distribution. 

The NFW halo is said to be universal, but when applied 

to rich galaxy clusters, it fails to explain the mass distribu- Astronomical Society Meeting, 200, 000 

tion; a nalysis of galaxy-galaxy correlations bv lSeldner fc Pee bles Bullock, J. S., Kolatt, T. S., Sigad, Y., SomerviUe, R. S., 
(1977) showed that the distribution of galaxies (n ga i <x A - v -> K1 yP in ' A - A - Primack > J - R > & Dekcl - A - 2001 

R~ 2A ) around rich clusters extends to 15 h _1 Mpc, with- 
out a feature that would suggest the expected sharp de- 
cline in density beyond R ma x (~ 0.8 Mpc), or at th e viria l 
radius (~ 2 Mpc). Recent results of iMcKav et all l)2002j) 
substantiate this claim using isolated groups of galaxies. 
They find n ga i cx r , and a LOS velocity dispersion er r 
that shows no significant variation between apertures of 



results appear to support an isothermal mass distribution 
- either truncated at a very large radius, or untruncated. 

Evolved ISOT baryon profiles whose halos were assumed 
to be un-evolving, were found to be more compact, and 
tended to form stars at velocities v c J> 25 kms -1 , leav- 
ing a gap between the largest model clouds and smallest 
galaxies and CHVCs. "Grown" halos, however, allowed 
clouds with circular velocities up to v c ~ 35 kms" 1 with- 
out forming stars, virtually closing the gap between model 
clouds and observed galaxies. The profiles which result 
from a gradually growing halo produce a significantly bet- 
ter match to the shape of the CDS, though the excess at 
low EW does, and should remain because the sub-halos 
are still there affecting the cloud cross section to absorp- 
tion at low column density (though their effects cannot be 
modeled with a 1-D simulation). 
I make the following conclusions: 

1. N-body/hydro simulations fail to predict the ob- 
served void clouds; diffuse sheets cannot explain the 
absorbers. 

2. The NFW halo must be rejected as an absorber 
model on a similar basis; it cannot explain the ob- 
served void clouds if "NFW" is taken to adequately 
represent the mass distribution beyond the virial ra- 
dius. 

3. A halo more massive than NFW is needed to explain 
observations. I find that a model consistent with a 
non-singular inverse square DM density profile, with 
a large truncation radius, is consistent with the data. 
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